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Abstract 

We prove the renormalizability of the generahzed Edwards model for self-avoiding polymer- 
ized membranes. This is done by use of a short distance multilocal operator product expansion, 
which extends the methods of local field theories to a large class of models with non-local singular 
interactions. This ensures the existence of scaling laws for crumpled self-avoiding membranes, 
and validates the direct renormalization method used for polymers and membranes. This also 
provides a framework for explicit perturbative calculations. We discuss hyperscaling relations 
for the configuration exponent and contact exponents. We finally consider membranes with long 
range interactions and at the 0-point. 
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1. Introduction 

Polymerized or tetliered membranes and tlieir statistical properties are a rich and 
interesting subject [|l],|^ , still under investigations. These objects can be viewed as a simple 
two dimensional generalization of one dimensional polymers Q . However such membranes 
exhibit a larger variety of behavior than polymers, due to the relevance of bending rigidity 
for two dimensional objects (as can be seen by simple dimensional analysis) . In 

particular they might undergo a crumpling transition, separating a high temperature (low 
rigidity) crumpled phase from a low temperature (high rigidity) flat phase with anomalous 
elasticity [|5|-|T0| . Also different behaviors are obtained if the membrane is fluid, i.e. without 
shear modulus. 

In this paper we consider only the case of tethered (polymerized) membranes, with 
internal in-plane elasticity, and we focus on the effect of self-avoidance constraints on the 



crumpled phase. For that purpose, a generalization of the celebrated Edwards model ^1 
for self-avoiding polymers has been introduced several years ago in [|T^,|T^ . This model 
incorporates the internal elastic properties of the membrane (in the crumpled phase) by a 
simple Gaussian entropic term, and self-avoidance through a local contact interaction. It 
allows for a systematic perturbative expansion in the steric interaction, and for a renormal- 
ization group treatment similar to the direct renormalization used for polymers within the 
Edwards model. At zero order, the crumpled membrane is described by a free Gaussian 
term, leading to a size exponent v given by 

- = — , (1.1) 

where D is the internal dimension of the object {D = 1 for polymers, D = 2 for mem- 
branes). Dimensional analysis shows that the canonical dimension for the steric interaction 
is 

2 — D 

e = 2D —d , (1.2) 

where d is the dimension of the bulk space in which the membrane is embedded {d — 3 for 
most physical situations). Taking e = sets the upper critical dimension dc 

AD 

dc{D) = ^—^ (1.3) 

below which self-avoidance is relevant (in the usual renormalization group sense). For 
polymers dc{l) = 4 and a systematic e = 4 — d expansion can be performed, using pertur- 
bative renormalization group techniques |1T^ - |T6[] . However for membranes dc{2) = +oo. 
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and such a perturbative expansion cannot be performed directly at D = 2. Still the 
model can be considered at a first stage for D-dimensional membranes with internal di- 
mension < D < 2, and renormalization calculations have been performed at first order 
T^ , p!3| , p!7| . [T^] , by expanding around the critical dimension both with respect to D and d. 



Thus starting from a D = Dq < 2 and d = dc{Do), one can reach the physical line D = 2, 
d < oo describing two dimensional membranes. 

The calculations within this model are parallel to those of the direct renormalization 
for polymers [13, TBI . The consistency of these calculations for polymers (in other words 



the renormalizability of the Edwards model) is ensured by the famous equivalence ^9 



between the Edwards model and a local quantum field theory (QFT) in the (i-dimensional 
embedding space: the 0(n) model in the limit n — 0, and by mapping the results of 
renormalization group theory for local QFT [^,21| onto the Edwards model [22,23] . Such 
an equivalence breaks down as soon as D ^ 1. In order to validate the previous calculations 
for membranes, an important issue is therefore the renormalizability of the generalized 
Edwards model for arbitrary D. This is not a simple problem. Indeed, although local 
in the embedding space, the interaction is non-local in the internal L>-dimensional space, 
since contacts may involve points arbitrarily far apart along the membrane. Thus the 
Edwards model can be viewed as a non-local field theory in internal space. For general 
non-local theories no theory of renormalization exists. 

A first step in this direction was made in [^, where it is shown that the theory is 
renormalizable at first order in perturbation. In this paper we establish the renormaliz- 
ability of this theory to all orders in perturbation theory. This is done by using previous 
results by the authors for a simpler model [^ , together with a new short-distance multi- 
local operator product expansion (MOPE), which extends methods of local field theories 
to a large class of models with non-local interactions, to which the generalized Edwards 
model belongs. A summary of these results has already been published in [^, and they 
are presented here in much greater details, with some applications. Our results validate the 
direct renormalization scheme used before. They also ensure the existence of renormaliza- 
tion group equations leading to scaling laws for self-avoiding membranes in the crumpled 
phase. They finally provide a general and practical framework for explicit calculations, 
which have been used in particular in where explicit O(e^) estimates for the size 

exponent u for self-avoiding membranes have been obtained, and in |2^ where the behavior 
of tricritical membranes at the 0-point is studied more thoroughly. 
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The paper is organized as follows. In Section 2 we recall the definition of the model 
and its perturbative expansion. This expansion admits a diagrammatic representation in 
terms of " dipoles" interacting via a Coulomb potential in internal L>-dimensional space. A 
proper definition of the model in non-integer dimension D is obtained by use of " distance 



geometry", following a previous work p5[. Finally, we address the problem of infra-red 
(I.R.) divergences, first showing that observables which are invariant under global trans- 
lations in embedding space are I.R. finite at zero order. 

In Section 3 we focus on short distances, or ultra-violet (U.V.), divergences, by first 
identifying those configurations which are singular, and by showing that they always cor- 
respond to a limit of short distances between some end-points of the dipoles. Around these 
configurations, we derive a multilocal operator product expansion (MOPE), involving gen- 
eral multilocal operators, which fully encodes the behavior of the amplitudes for the model 
at short distances. This MOPE is then used to classify the U.V. singular configurations 
by power counting. We show that these U.V. divergences are proportional to the insertion 
of multi-local operators. At the upper critical dimension, i.e. at e = 0, the only dangerous 
operators are those present in the original Edwards model, which strongly suggests that 
the model is renormalized onto itself. 

Section 4 is devoted to the proof of renormalizability of the model: U.V. divergences 
can be subtracted and the theory can be made U.V. finite at e = by adding appropriate 
counterterms to the Hamiltonian of the Edwards model. We focus on the simple case of 
infinite membranes with internal flat geometry. We first show that the model is U.V. finite 
for e > (super-renormalizable case) (Section 4.1), as well as I.R. finite when e — > 0-(-, as 
long as translationally invariant observables are concerned (Section 4.2). We then analyse 
the U.V. divergences at e = 0, making repeated use of the formalism previously developed 
in for the model of a free membrane interacting with a single point 13^,0 . Indeed, 



this formalism, developed for local but singular interactions, can be adapted to the class 
of singular multilocal interactions that we consider here. In Section 4.3, single U.V. singu- 
larities (superficial divergences) are classified in terms of dipole configurations which are 
called diagrams, following [^. A subtraction prescription for these superficial divergences 
is then introduced. In Section 4.4 the nested structure of multiple U.V. divergences and 
the required subtractions are presented on simple examples at second order. Iterating this 
analysis naturally leads to a formulation of renormalization in terms of nested families of 
successive divergent subdiagrams, which are analogous to Zimmermann's forests of renor- 
malization theory in perturbative local QFT. A forest formula for a subtraction operator 
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is given in Section 4.5, and is shown to correspond to the introduction of multilocal coun- 
terterms in the Hamiltonian of the model. Section 4.6 is the core of the proof. We show 
that the renormahzed observables are U.V. finite at e = 0, by using a strategy analogous 



to that of The distance integration domain is decomposed in sectors, analogous to 
Hepp's sectors in perturbative local QFT, and forests of subtractions are reorganized with 
respect to each sector in order to ensure U.V. finiteness. The argument is quite techni- 



cal, and relies on the techniques developed in [g5[, but is presented here in a consistent 
way. Finally in Section 4.7 we discuss some additional points, in particular we show that 
renormalization does not affect the I.R. finiteness of invariant observables at e = 0. 

The rest of the paper presents applications of this renormalization formalism. In 
Section 5 we derive the renormalization group equations for the model, and we show how 
the scaling laws for the self-avoiding membrane model are derived for e > 0. As an 
application, we perform renormalization at first order in perturbation theory, using our 



formalism, and we show how the 0{e) result of |T2| , p!3[1 for the size exponent ly is recovered. 

Section 6 is devoted to the important case of finite membranes. Indeed, the direct 
renormalization method requires the study of finite membranes (or polymers), in partic- 
ular since it is the internal size L of the membrane which fixes the renormalization scale. 
In Sections 6.1-4 we show in details how our formalism can be extended and applied to 
the case of finite (open or closed) membranes, possibly with non-fiat internal geometry 
(i.e. with frozen Gaussian curvature), and to the case of several interacting membranes. 
We show in particular in Section 6.5 how to derive finite size scaling laws for finite mem- 



branes, and we discuss the hyperscaling relations [T^] between the size exponent u and 
the configuration exponent 7. In Section 6.6 we discuss in details the relation between 
the direct renormalization method and our renormalization formalism, and we prove that 
the scaling assumptions which are the starting point of direct renormalization are indeed 
correct. Finally in Section 6.7 we discuss the contact exponents which describe the scaling 
laws for contacts between different elements of the membrane. In particular we show how 
to take into account properly the edge effects in the calculation of edge contact exponents. 

Finally Section 7 is devoted to two other models of interacting membranes, which are 
also amenable to our formalism. The first case is fiexible tethered membranes with long 
range repulsive interactions, for instance uniformly charged membrane with non-screened 
Coulomb interactions. We show that for long-range interactions our formalism implies 
that there is a renormalization of the internal elasticity, but no renormalization of the 
interaction (charge), and that this leads to the fact (already well known for polymers) that 
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the size exponent v is exactly that given by a Gaussian variational method. The second 
case is that of membranes with competing attractive and repulsive short range interactions, 
corresponding to the generalization of polymers at the G-point (separating the swollen 
phase from the collapsed phase) . While for polymers it is well know that a repulsive three- 
body interaction is relevant to describe the scaling properties of the 0-point (at least in 
the £ = 4 — (i expansion) , we show that for membranes there is a competition between this 
three-body interaction and a modified two-body interaction. We then perform analytically 
the first order calculation of the renormalization group fiow for this new interaction. We 
refer to ||2^ for a subsequent and more complete analysis of this problem and a treatment 
of the crossover between the two interactions. 

Section 8 is the conclusion and discusses still open problems. Some more technical 
points are discussed in the appendices. 



2. The Perturbative Expansion 

2.1. Perturbation Theory Diagrammatics 

In this section we define the model by giving the formal rules of its perturbative 
expansion. We start from the Edwards Hamiltonian, generalized to the case of a manifold 
with arbitrary internal dimension D. Each point on the manifold is labeled by a D- 
dimensional vector x G IR^ . The position of the points x in the embedding d-dimensional 
Euclidean space defines a vector field r(x) G IR*^. The continuum Hamiltonian reads 

^[il/keT = ^ Jd'^x (V,r(x))' + ^Jd''xJ d^x' 5''{r{x) - r{x')) . (2.1) 

The first Gaussian term (Vcr(x))^ = ^^^i{dr{x) / dx"')'^ describes the local elastic energy. 
This term alone describes for 6 = a "phantom" polymerized membrane with no self- 
avoidance, and is mostly of entropic origin. The second term is a two-body short range 
5-potential which describes weak self-avoidance with excluded volume parameter 6 > 0. 
This term is both singular (it is a singular function - in fact the Dirac distribution in IR*^ 
- of the field r) and non-local in the internal space IR^ . In this respect this theory is quite 
different from usual local field theories. With our choice of units, the Hamiltonian (2T) is 



^ In this paper, bold quantities like r = {r'^}, fi = 1, . . . ,d refer to the external d- dimensional 
space, and non-bold quantities like x = {x"}, a = 1, . . . , D to the internal D-dimensional space. 
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dimensionless. The Hamiltonian (2J.) is invariant under global translations in embedding 
space r(x) r(x) +ro. In ( |2.1D the integrals run over IR^, which means that we consider 
an infinite membrane. The case of finite membranes will be discussed in Sec. 6. 

The partition function is defined as a sum over all membrane configurations r(x). 
More precisely, in this continuum formulation it is given by the functional integral: 

Zt = Jv[r{x)]exp{-H[r]/]^^T) . (2.2) 

It has a perturbative expansion in b, formally given by expanding the exponential of the 
contact interaction 

ife) iV! J .=1 \=i 

where Zq is the partition function of the Gaussian manifold, and (■ • ^ denotes the average 
with respect to the Gaussian manifold (6 = 0): 

((■ ■ ■))o = ^ / [-\ I d^'x {V^r{x)f^ (•■■). (2.4) 

Physical observables are provided by average values of operators, which must be in- 
variant under global translations. Using Fourier representation, local operators can always 
be generated by the exponential operators (or vertex operators), of the form 

V^(z) = e^^-'^^^) . (2.5) 

In perturbation theory the field r(x) will be treated as a massless free field and the momenta 
q will appear as the "charges" associated with the translations in IR^. Translationally 
invariant operators are then provided by "neutral" products of such local operators, 

p p 

1=1 1=1 
The perturbative expansion for these observables is simply 

(n^e^qrr(.,)) = — ^ ^ ^, ^ y rf^X, ( e^^' •'■(^') 5'^(r(x2a) - ?(x2a-l )) >o • 

(2.7) 



+k 
• — 



-k 



Hk. 1 




fig. 2 

Each 5 function in (|2.3| ) and ( |2.7| ) can itself be written in terms of two exponential operators 
as 



5%r{x2)-r{x^)) 



5^(ki +k2)e^''^-^(^i)e'''2-^(^2) 



2.^ 



Viewing again the momenta ki, k2 as charges assigned to the points xi, the bi-local 
operator ( |27^ ) corresponds to a dipole, with charges ki = k, k2 = — k, integrated over its 
internal charge k. We depict graphically each such dipole as on fig. |l|. 

Similarly, the product of bi-local operators in ( p.3| ) and (|2.7| ) can be written as an 
ensemble of dipoles, that is as the product of 2N vertex operators with "dipolar 
constraints" 

CjkJ = (27r)^5'^(k2a-l+k2a) , (2.9) 



then integrated over all internal charges k^: 



N 



2N d'^ki N 



2N -r 



n 5''{rix2a)-rix2a-i)) = I n n Ca{k,} U e^'^^''^^^^ 

a=l J i=l [2itP a=l i=l 



(2.10) 



Products of such bi-local operators and of external vertex operators, as in (|2.7|) , are de- 
picted by diagrams such as that of fig. ^ 



The Gaussian average in (|2.3|) is easily performed. We use the identity 
<ne^^^-(-)>^ = exp(-i^k,.k,(r(x,)r(x,))o 



(2.11) 



where r is any single component of r. Thanks to the neutrality condition k^ = 0, we 
can rewrite it as 



1 



^j^gik..r(x.)>^ = exp(--5^k,-k,G'(x,-x 



(2.12) 



with the translationally invariant two-point function 

1 2 
G{x^,Xj) = -- ((r(x,) -r(xj)) )o 



|2-D 



i2-D)SD 



Here 



S 



2 7r^/2 



D 



(2.13) 



(2.14) 



r(I^/2) 

is the volume of the unit sphere in D dimensions. In ( |2.13| ) and ( |2.14| ) the standard analytic 
continuation in the dimension D to < -D < 2 is performed. 

Integration over the momenta k^ then gives for the A^'th term of the perturbative 



expansion for the partition function Zb (|2.3|) the "manifold integral" 

2N 



{2n) 



-Nd/2 



n d^XiA{xi}~^ 

i=l 



[2.15) 



with Ajxi} the determinant associated with the quadratic form (on IR ) Qjki} = 

2N 

S kikjG'(a;i, Xj) restricted to the A^-dimensional vector space defined by the neutral- 

ity constraints Ca{ki}, k2a + k2a-i = 0. A{xi} is given explicitly by the determinant of the 
N X N matrix Aab (with row and columns labeled by the dipoles indices a,b = 1, ■ ■ ■ , A^) 

A = det (Aafc) , Aab = G{x2a-l,X2b-l)+G{x2a,X2b)-G{x2a-l,X2b)-G{x2a,X2b-l) 

(2.16) 

Similarly, the A^'th term in the perturbative expansion of the P-point observable ( |2.7| ) 

is 



(27r) 



-Nd/2 



2N 



n d^Xi A{xi} 2 exp ( - ^ E ■ q, 



i=l 



1 P 



2 Lm=l 



N^"^ is the (/m) minor of the (P + A^) x (P + A^) matrix 



G{zi,Zm) G{zi,X2b-l) - G{zi,X2b) 

G{X2 ) - G{X2 



l<l,m<P 
l<a,b<N 



(2.17) 



(2.18) 
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2.2. The Distance Measure 



In the integrals ( |2.15| ) and ( |2.17| ) the dimension d of external space appears only as 
a power, and can be considered as a continuous parameter i < (i < oo. For the reasons 
explained above, we also need to treat the internal dimension D as a continuous variable 
< D < 2. This can be done, thanks to the formalism developed in [^. The integrand 
in (|2.15 ) is a well defined function of the relative Euclidean distances between the internal 



points dij = \xi — Xj\, thanks to the Euclidean invariance of the model in internal D- 
dimensional space. The integral over the M = 2N points Xi in IR"^ can be rewritten as an 
integral over the squared distances 



Gij = dij 



[Xi Xj^ 



(2.19) 



times a global internal translation factor V = J-^d d xi (the volume of the membrane): 

M=2N 



n 



D 



Xi 



•^M l<i<j<M 



a 



(2.20) 



with the measure on the integration domain Am (defined below) 



_ M(M-l) 

2 2 Sd Sd-1 



Sd-m+2 

where PM[aij] is the Cayley-Menger determinant 



{D-M)/2 



(2.21) 



(-1) 



M 



}M-1 






1 


1 


1 


1 





ai2 ■ 


■ ■ aiM 


1 


a2i 





■ 0,2M 


1 


ttMl 


aM2 ■ 






(2.22) 



The integration domain Am is the subset of the variables a^j which can be realized as 



actual squared distances between M points in IR^''^"^. It is given by the M — 1 constraints 
> 0, 2 < i^T < M i. Geometrically the Cayley-Menger determinant is simply the 
squared volume of the (M — l)-dimensional parallelepiped generated by the vertices xi, . . ., 



^ Strictly speaking the observables ( |2.17| ) are then defined as functions of the invariants -qm- 
For instance in the case of three points (M = 3), the two non-trivial constraints P2 > 
and P3 > imply that the squared distance variables aij are positive, and satisfy the triangular 
inequahties | Ja\2 — -y/ois | < a/ois < 1/012 -|- 1/023 • 
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xm- Originally the correspondence ( 2.20 ) is valid only for integer dimension D > M — 1. It 
can then be analytically continued for non-integer dimension < D < oo. For D > M — 2 
the measure ( p.21|) remains integrable. For D < M — 2, (|2.21j ) becomes singular due to the 



vanishing of the sphere volumes Sd-k at D = iiT — L, L G IN, i.e. for integer D < M — 2, 
and to the divergences of Pm[ci]~^^~'^'^^^'^ when the determinant PM[fi] vanishes. It has 



been shown in that these divergences occur on the boundary of Am, and can be dealt 
with for non-integer D hj a, finite part integration prescription. Thus /U^[a] becomes for 
D < M — 2 a well-defined distribution with support on Am- For integer D < M — 2 
this distribution concentrates on the boundary of Am, and more specifically on those 
submanifolds which describes distances which can be realized as actual distances between 
M points in R^. 



For correlation functions like ( |2.7| ) one has to perform the integral ( p.l7|) . The inte- 



grand is now a function of the squared distances aij, an and aim between the M = 2N 
internal points Xi and the P external points zi. The integral over the internal points is 
rewritten as an integral over the variables a^j and an, the distances between the external 
points, that is the aim, being kept fixed: 

M 

/ JJfi^Xi = / Yl ^^ij Yl dan fx^p[aij,aii,aim] , (2.23) 

"^^^ i=l J Am,P i<i<j<M l<i<M 

1<1<P 

with the measure depending now on the Cayley-Menger determinant Pp[aim] relative to 
the P external points, and on the Cayley-Menger determinant PM+p[aij, an, aim] relative 
to the M + P external and internal points: 



jy M{M+2P-1) \ ^ IVl-tf I'^lJ, ^IL, ^Im 

fJ'M,p['^ij:'^ih'^lm] — 2 2 S D- pj^i . . . S D- p- M +2 



(^PM+p[aij,aii, ai 

r 



D~M-P 



D-P 

\ 

ml 



Pp[air 

(2.24) 

Am,p = AM,p{aim) is now the subset of the variables {aij.an) which can be realized as 
actual distances between M points and between these M points and the P fixed external 
points in IR^"*"^"^. For generic aim va. Ap, this measure is integrable for D > M + P — 1, 
and becomes a distribution for D < M + P — 1. Alternative parametrizations of the 
measures ( |2.21| ) and ( |2.24| ) are given in ||25|| . 
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2.3. Infra-Red Problems 

The perturbative expansion of the model has potential infra-red divergences for D < 2. 
In fact, infra-red divergences arise already for the free theory (6 = 0) if we deal with observ- 
ables which are not translationally invariant. Indeed, the free massless propagator which 
appears in ( |2.11| ) is infra-red singular for D <2. It can be defined properly by introducing 
an infra-red regulator (represented by a length scale L) , for instance by considering a finite 
membrane @ with internal linear size of order L, or by adding a mass term J m?r^ in the 
Hamiltonian with m oc 1/L, corresponding to a confining potential in external space. 
With such regulators, the free propagator is finite and behaves, for L | 

For translationally invariant observables, we have seen that the neutrality conditions over 
the external momenta and the neutrality constraints Ca imply that in ( |2.11|) the infra-red 
divergent propagator can be replaced by the translationally invariant two-point function 
G{xi,Xj) given by ( |2.1tj| ), which is infra-red finite. In fact G{xi,Xj) can be considered 



as the infra-red finite part of the free propagator ( |2.25| ), defined for instance by letting 
L (X) for D > 2, and then by analytic continuation to D < 2. This is a first example of 
cancellation of infra-red divergences. 

Another source of infra-red divergences comes from the integration over internal points 
Xi in ( p.l5|) and (|2.17|) , when some of these points go to infinity. However, we expect these 



divergences to disappear when considering average values like (|2.7|) , by a compensation 
of these divergences in the numerator, given by integrals like ( |2.17|) , and those in the 
denominator Zf,, given by integrals like ( |2.15| ). At variance with the cancellations of I.R. 
divergences in the integrand, which occur for < -D < 2 and for any < d < oo, these 
cancellations should occur only in the vicinity of the upper critical dimension dc{D) = 
4D/{2 — D) where the model is renormalizable. For instance, one can check by an explicit 
calculation that at first order in 6, the invariant observables ( |2.7|) are I.R. finite as long as 
e = 2D — (2 — D)d/2 < D. Similarly, the free energy density f = — \n{Zi,)/V (V being the 
volume of the membrane) is I.R. finite at order b, provided that e < D. At this stage, we 
can only make the following conjecture: 



This requires a proper treatment of the global zero mode, see 
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The term of order b'^ of the perturbative expansion of invariant observables is I.R. 
finite for < D < 2 and e < D/n. 

This statement is analogous to the Ehtzur theorem p3| , |3^ , which states that invariant 
observables are IR finite in two dimensional non-linear sigma models. We shall come back 
to this point later. 



3. Ultra- Violet Singularities and the Multilocal Operator Product Expansion 

3.1. UV Singular Configurations 

A first condition for the manifold integrals ( |2.15|) and ( |2.17| ) to be well defined is that 
the determinant Ajx^} must be nonnegative. This can be shown by methods used to prove 



the so-called Schonberg theorem |3^, that we now discuss, using a convenient electrostatic 
formulation. 

Let us consider the quadratic form 

M 

Q{kJ = k^k,G'(x„x,) (3.1) 

restricted to the (M — 1) -dimensional subspace E C defined by the global neutrality 
constraint 

M 

J]k, = . (3.2) 

i=l 

Here G{xi,Xj) is given by ( p.l3|) , so that (5{ki} can be viewed as the electrostatic energy 
of a globally neutral system of M charges k^ in IR^. We still take < D < 2, and make 
use of distance geometry to view Q as a function of the squared distances a^j = (x^ — Xj)^, 
taken in the domain Am^ which can be realized as actual squared distances between M 
points Xi in IR^ , with M' > M — 1 integer. 

Schonberg Theorem: 

For < D < 2, the quadratic form (5{kj} is positive on E. Moreover, (5{ki} = if and 
only if the local "charge density" (defined in IR^ , M' > M— 1, by p(x) = J^fLi ki5(x— x^)) 
is zero everywhere. 

Proof 
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We simply use the Fourier transform representation of the kernel —\xi — Xj\'^'^ in IR 
valid for < < 1, 



121^ 



C{M\iy) 



(27r)^ 

with C{M', v) a positive constant (for < u < 1 



7 m 



2u-M' ^^iq{xi-Xj) _ 



(3.3) 



C{M',u) 



M' ■ 



TT 



Taking = (2 — D)/2 and using the neutrality constraint = we get: 



(3.4) 



g{kj 



C{M',v) 

(2 - D)Sd J (27r) 



_^ ^ \ \-2u-M' 



M 



ki e^'?^' 



(3.5) 



Since the kernel \q\ '^^ ^ is positive, the positivity of Q follows. Moreover, Q = if and 
only if the Fourier transform of the charge density p{x) 

M 

p{q) = ^kie^^^» (3.6) 



1=1 



is zero, i.e. if p vanishes everywhere. 

We now consider the quadratic form Q associated to the M = 2N points Xi in ( |2.15|) . 
Since Q{ki} is positive on the (M — l)-dimensional subspace E, it is also positive on the 
A'"-dimensional subspace F C E defined by the dipole neutrality constraints Ca{k-i} 



F = {{kaclR'^: k2„-i-k2„ = 0, a = l,...iV} 
This implies that A{xi} = detpiQ) > in A2n- Similarly, the quadratic form 



(3.7) 



(3.^ 



l,m=l 



in the P-point integral ( |2.17| ), which depends on the P external charges qz, such that 
^^j^ = 0, is given by the minimum over all possible internal charges {k^} G -F of the 
global quadratic form Q{qi,ki}. Therefore Qiqi} > 0. 

Still the integrals ( |2.15| ) and ( p.l7|) may diverge if A{xi} = 0, i.e. if and only if the 
quadratic form Q has a non-zero isotropic vector {k^} G F, such that 



Q{ka = , {ka^{o,...,o} . 



(3.9) 
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fig. 3 

The Schonberg theorem aUows us to identify precisely when this happens: A{xi} = if 
and only if the 2N dipoles have their positions such that the charge density p{x) can be 
kept zero everywhere, while some dipoles have (appropriate) non-zero charges. If all the 
points i occupy distinct positions in IR^, this is clearly impossible. Thus Ajx^} = first 
requires that some points must coincide. More precisely: 

(1) To any singular configuration is associated a partition of the 2N points into "atoms" 
V, such that the points i which belong to the same atom V share the same position 
x-p (in the distance space this means i,jEV^ aij = 0). Some atoms V may consist 
of only one point {i}. Since the points i are end-points of the dipoles a, the atoms 
V are assembled into "molecules", defined as connected sets of atoms attached by 
dipoles (see fig. A singular configuration is thus associated to a grouping of the 
points into molecules. 

Moreover, the condition of zero charge density means that each atom V must be neutral: 

^k, = , (3.10) 

while at least one dipole carries a non-zero charge. Clearly this is possible if and only if 
the non-zero dipoles belong to a closed loop in the molecule. Thus A{xi} = requires in 
addition that: 

(2) One of the molecules contains at least one closed loop of dipoles. 

In fact, for a given molecule, the non-zero charges k^ are carried by the one-particle irre- 
ducible (IPI) components of the molecule, that is by the connected subsets which remain 
connected when any line (dipole) is removed. 
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3.2. The Multilocal Operator Product Expansion (MOPE) 

We now analyse the short distance singularity associated with a given singular con- 
figuration, i.e. a molecule characterized by a set of atoms V with positions x-p, and by a 
set of dipolar links a between these atoms. 

For this purpose, we consider the points i which belong to a given atom V and see how 
the integrands A~'^/^ and A~'^/^e~^'^ behave when the points i collapse inside V {xi — > xp, 
i E V). The singularity can be analyzed in distance space, by letting the squared distances 
aij go to zero for i and j belonging to the same atoms. A more convenient and physical 
way to study this behavior is to formulate the problem in terms of the "physical" positions 
Xi of the points i and x-p of the atoms V. Let 

Hi = Xi — x-p , for i E V (3-11) 

be the relative position of point i with respect to its atom V. 

The short distance singularity of A"*^/^ is analyzed by performing a small yi expansion 
of the product of bilocal operators 

n 5\v{x2a)-r{x2a-i)) = / H H CJkJ n e''^-'-^^^) (3.12) 

a=l J i=l (iTTj" a=l i=l 

in the Gaussian theory. This is done as follows: 

(1) We regroup the vertex operators in each atom V. For each atom we factorize the 
product of vertex operators into its expectation value times its "normal product" 

n = (n e^^»-'(^»))o: n e^'^^-''^^^): . (3.13) 

iev iev iev 

This can be considered as the definition of the normal product : :. At the diagram- 
matic level, this means that when evaluating the expectation value of the l.h.s. of 
( ^.13| ) times the other "external" operators, using Wick's theorem (which leads to 
( ^.11| )) one can always factorize out the contribution involving pairs of points («, j) in 
the atom V, which is nothing but the average value: 

e '-^^^ = ( n e'''-''("^))o . (3.14) 

iev 

What remains is the contribution involving propagators between a point i in the atom 
V and the other points not in V. This is precisely the definition of the normal product, 
which consists in omitting internal contractions between points inside the : :. 
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(2) Now the normal product relative to the atom P is a smooth function of the positions 
Xi in a neighborhood of x-p, since the short distance singularities at Xi = Xj have been 
factorized out. Therefore we can perform a Taylor expansion of this normal product 
in the yi, whose coefficients are normal ordered local operators located at x-p, and 
which is of the form 



-Q giki-r(xi). _ n e^' n e^'^*'''*-^^'' 
iev iev iev 



Xi=X'p 



^ . ^ i( E k,)-?(^p) (3-15) 

= J2C''{y^.K}■.A{xJ,)e 

A 



The local operators A (x-p) are monomials in the derivatives of the field r at the 
point x-p. The corresponding coefficients C^{yijhi} are monomials in the yi and 
variables, for i E V. For instance, the first terms of this Taylor expansion are given 
explicitly by: 

: n e'^^'"'^'''^: = :e ^^t^ : + i S y^kf :Vx,r''(xp)e : 

iev iev 

- - E y,2/,kfk^^:V,,r^(xp)V.^.r'^(xp)e : (3.16) 

^ ijev 

i „ i( S ki)-?(a;-p) 

+ - S ?/i?/ikf :V,y,,r'^(xp)e : + ■■■ 

/ iev 



where the summation over the D components of the y^'s are implicit: y^Vx 



(3) From translation invariance in internal space IR^, the average value ( ^.14| ) relative 
to the points z G P is the same when evaluated in terms of relative positions, i.e. 
{r{y,)r{yj))o. 

(4) Thanks to the global neutrality condition on the k^ in (|3.12|) , we can replace all 



propagators by their infra-red finite part ( |2.13| ) . This means on the one hand that the 
correlation function (r(j/i)r(j/j))o of point (3) above is taken to be G{yi, yj), building 
the factor 

-i E ii-iiG{yi,yj) 

e '■^^'^ (3.17) 



in ( p.l4| ). On the other hand, the operators in ( |3.15| ) are also defined through the 



normal order prescription, which uses the I.R. finite propagator ( |2.13| ). This means 
that when evaluating correlation functions of products of such operators, located at 
positions xp, xp', . . . one uses the propagators G{x'p, xp') and their derivatives. 
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(5) We now gather the contributions relative to the atoms V and hnks a of a given con- 
nected molecule M.. If the singular configuration consists in more than one molecule, 
the following treatment can be applied independently for each molecule M.. For each 
atom V we introduce its total charge k-p = ^j^pk^. This is done by inserting 
1 = J (i^k-p 5'^(k-p — Eie-pki) in (|3.12|) , which becomes an integral over the k^ and 
the kp>. The constraints Ca for the dipolar links a in the component M. ensure that 
YliVeM ~ ^- therefore rewrite 

n C„{kJ = n' C„{kJ(27r)'^5'^( E kp), (3.18) 

aeM aeM veM 

where II' means that one of the constraints has been omitted arbitrarily (the result 
does not depend on the link a which is omitted). 

(6) We now Taylor expand each 5'^(kp — Eiep k^) in powers of kp: 

5^(kp- E k.) = Y^^^^i^^t^'^- S k.) (3.19) 

m 

where the sum runs of (i-uples m of non-negative integers m = ... d, with the 

standard notation: 

A n-.ni,. ^m ^ 9 \ ""^ ^ (mi + --- + md)! 

m! = 



mi! ■ • • m^! 

(3.20) 

(7) At this stage the product of bilocal operators ( |3.12| ) has been expanded as a sum of 
terms of the form 

where 

rf^k,- 



c{A.,r^.}|^j = / n ^ n; cjkjx 

J ieM [Zn)'^ aeM 



X n <^C^^{2/„kJ,eP^(Vj^)'"^5^(- E kOe -^^^ 
VeM I m-p! iev 

(3.22) 

and 

^{A^,A„}{xv]^ I n d^^v {'i'^Y 5\ kp) n |(kp)--:Ap(xp)e^'^-''(^-) 
J VeM VeM VeM I 

(3.23) 
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The sum extends over all collections of {A-p^m-p} relative to the atoms V of the 
molecule A^. Once integrated over the k^, (|3.22|) builds a function of the relative 
positions yi inside the atoms. Since G{yi, yj) is a pure power of \yi — 2/j |, this function 
is moreover homogeneous with respect to the yi . The integral over the h-p in ( |3.23| ) can 
be performed by performing a Fourier transform of the 5-function, thus introducing a 
variable r: 



${A^,m^}{xp} = yrf'^r ^n^|:Ap(xp)(iVp)'^^5^(r-r(xp)):| . (3. 



24) 



Eq. (|3.24|) describes generically a multilocal operator, involving contact interaction (at 
point r) in external space between the points x-p, associated with the different atoms 
V of the molecule M., with local operators Ap{xp>) at each point x-p. As can be seen 
in ( |3.15|) , these local operators, which must be translationally invariant, involve only 
derivatives of r(x). 

Multilocal Operator Product Expansion (MOPE): 

The result of points (l)-(7) above is the following multilocal operator product ex- 
pansion property: When approaching the singular configuration characterized by the con- 
nected molecule A4, the product of bilocal operators corresponding to the dipoles a of the 
molecule can be expanded as a series of multilocal operators $ 

n (5^(?(x2a)-?(x2„-i)) = V C{^^'*^>{y,} ■ (3.25) 

The $'s are multilocal |A^|-body contact operators, where is the number of atoms 
in A4. They depend only on the positions xp of these atoms V. The coefficients C of 
the MOPE depend on full structure of the molecule A4, and on the relative coordinates 
{Ui = Xi — x-p}i^-p of the end-points inside each atom V of separately. They do not 
depend on the relative coordinates {xi — Xj} of end-points belonging to different atoms 
{i E V, i E V, V 7^ V). Finally the coefficients C are homogeneous functions of the 
i/i's, with a degree of homogeneity deg[C] which depends on the molecule M. and on the 
{Av,Ap}. 
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A direct application of this MOPE is the short distance expansion for the integrands 
A~'^/^ and A~'^/^e~i'^ in ( |2.15| ) and ( |2.17| ). For instance, we have 



a=l 



{A-p ,m-p } ^ 



(3.26) 

where the expectation value of a product of bilocal operators times a multilocal operator 
in the r.h.s of ( p.26|) can be evaluated by the same rules as in Section 2.1, using the 



Fourier representation of (|3.23D for $. It is a function of the positions {xi}i^M ci-nd of the 



{x-p}'p^M- Using the Schonberg theorem, one can show that, like A this function is 
well defined and non-singular if the {xi}i^M ci-nd the {x'p}'p^M do not coincide. 

3.3. Power- Counting 

We are now in a position to analyze the singular configurations and the associated 
divergences. Given a connected molecule M., we contract the points i towards the center 
of the atoms V by rescaling the yt = Xi — x-p by a global factor A. The behavior of the 
integrand in (|2J5|) and (|2T7|) when A ^ is given by the MOPE. Since C^^^'^'^^lyJ 



is a homogeneous function of the {yi}, each term of the MOPE gives a single power of A. 
Prom (J3. 221) and using the fact that k scales as y~^, the degree of C'f"^''''"^^{yi} is given 



by 

deg[C{^^''^^>] = c^zy- (# link. -1) -rfzy-# point. + E {\inp\u - dim[Ap] + du) (3.27) 



aGA4 iGA4 atoms 



with, using (|3.2UD 



|m| ^ . (3.28) 

To obtain the corresponding degree of convergence uj, we have to integrate over the positions 
for all the points i but one inside each atom V, since the position x-p of each atom is fixed 
in the contraction process. This degree of convergence is thus given by 

^{A^.riip} _ deg[C^^^''^^>] +!)(# points -# atom.) (3.29) 

Using the fact that ^points = 2 #iinks and defining the number of internal loops of the 
connected molecule as 

#loops = #links — #atoms + 1 (3.30) 
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one obtains 

^i?""""^ = ^(#atoms-2)+e(#ioops)+ E ( |mp | z/ - dim[^^] ) (3.31) 

atoms 

where we recall that e and v are given by (|1.2|) and ( |1 . 1| ) respectively. 

The last sum in ( p.31|) , relative to the atoms "P, is always positive, since the local 



operators A-p generated in the MOPE are monomials in (multiple) x-derivatives of r, and 
have negative dimension (for L> < 2). Indeed, 

dim[^] = -#x-derivatives + Z^#r < (z/-l)#r = -y#r < (3.32) 

This sum is nothing but the dimension of the "dressing" by {A-p-, mp} of the molecule M.. 
The smallest degree is obtained by setting all the A-p = 1 (the identity operator) and 
all the mp = 0. This defines the superficial degree of convergence of the molecule M.: 

U)M = D (#atoms — 2) + e (#loops) (3.33) 

Only the molecules M. such that 

UJM < Q (3.34) 
give short-distance U.V. divergences in the integrals (|2.15| ) and ( p.l7|) . From the MOPE, 



we expect that the divergences given by a coefficient C^^'p^'^'p^ will be proportional to the 
insertion of the corresponding multi-local operator ^{A'p,rn-p}{'^'p}- Moreover, it is easy 
to check that the degree of convergence uj^^'"^^^ is related to the canonical dimension of 
the multi-local operator ^{A-p,m-p}{^'p}j integrated over the whole space, by 



+ ejj^ links (3.35) 



T^atoms 



as expected from dimensional analysis. 

Depending on the sign of e there are three cases: 

(1) e > 0: only the molecules such that #atoms = 1 and with a number of loops #ioops 
small enough diverge (#ioops < D(2 — #atoms)/e)- There is a finite number of such 
divergent molecules, and the model will be super-renormalizahle. 

(2) e = 0: only the molecules such that #atoms = 1 or 2, but with #ioops arbitrary, 
are divergent. There is an infinite number of divergent molecules, but we will show 
that they correspond to a finite number of multi-local operators. The model will be 
strictly-renormalizahle. Such divergent configurations are depicted on fig. || . 

(3) e < 0: all the molecules with a large enough #ioops diverge. There is an infinite 
number of divergent molecules, which corresponds to an infinite number of multi-local 
operators (in fact all of them). The theory will be non-renormalizahle. 
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(a) (b) 
fig. 4 



4. Renormalization 



We are now in a position to analyse the UV divergences of the theory at e = and to 
show that the theory is renormahzable, that is that the UV divergences can be subtracted 
and that the theory can be made finite, in perturbation theory, by adding appropriate 
counterterms to the Hamiltonian. We shaU foUow the general strategy of the so called 



dimensional regularization scheme, as used in local field theories and for the model of |25 
The theory is first defined for e > and small, and e is considered gularization 
parameter. Then we construct the renormalized theory by adequate subtractions and by 
taking the limit e ^ 0. 

4.1. The theory for e > 

When e > 0, the only divergent configurations are the molecules A4 made of one atom 
and of dipoles, with < D/e. According to the results of the previous section, the 
corresponding singularities are given by the MOPE coefficients C"^ when the dipoles 
are contracted toward a single local operator A{x) (we can here assimilate A{x) with the 
"1-local operator" m^oyi-^} defined by ( ^.241 ), and "1-local operators" with m 7^ are 
automatically equal to zero). The corresponding degree of convergence of each term 
is given by ( p.31|) , and is 

-D + eN-dim[A] (4.1) 



,A 
M 



We see that only the identity operator Aq = 1 with dim[l] = gives a divergence. Less 
relevant operators give > 0, for instance the least irrelevant one is (Vr)^ which has 
dimension dim[(Vr)^] = D, so that uj = eN > 0. 

Since for e > there is a finite number of divergent molecules, and since the associated 
divergences in the manifold integrals are proportional to the unity operator, this suggests 
that: 
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(1) Those divergences are cancelled by adding to the bare Hamiltonian ( |2.1| ) a counterterm 
A7i/kBT independent of r and of the form 

AH/kBT = A(6;e) j d^xl (4.2) 

(2) As a consequence, these short distance divergences should be absent from the manifold 
integrals when one computes in the interacting theory expectation values of physical 
observables O as defined for instance by ( p.6|) . Indeed the divergences, being inde- 
pendent of the configuration field r, factorize out of the functional integral J'E'ir] 
and should cancel between the numerator in ( |2.7| ) and the normalization denominator 
factor Zi,. This means that perturbation theory for physical observables in the infinite 
membrane model is ultra-violet finite as long as e > 0. 

Let us give the outline of the proof. For simplicity we denote by (j){x, y} the two-point 
operator 

Hx^y} = d''{r{x)-r{y)) . (4.3) 

and we consider the A^-th term of the perturbative series expansion of the expectation 
value of an observable 0{^.j{{zi}), as defined by (pT6|). This term is proportional to the 
integral of the connected expectation value (in the free theory) of O times the product of 
(p operators 

0^^^{z^} = / IT rf^x,(0{^j({z,}) n (^{X2a-i,a;2a})r" (4-4) 

The subscript {■ ■ ■)o refers to the expectation value in the free theory with 6 = 0. The 
connected expectation value (e. v.) (■ • of a product of A^ operators is defined here 

by inverting the well known relation 



{Ai---An) = T n ( n A.r^^^ (4.5) 

■i^subset 

partitions subsets of 

of { 1 , ■ ■ ■ , N" } the partition 

which gives 

(Ai-.-A^)^""^"^ = V ^(partition) TT { U A,) (4.6) 

partitions subsets of 

of { 1 , ■ ■ ■ , iV } the partition 

with 

(partition) = (-1)Ip^'^*I-i (|part| - 1)! (4.7) 
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where |part| is the number of distinct subsets (elements) of the partition. 

With these notations, in ( [4.4| ) the observable O has to be considered as a single 
multilocal operator, although it is a product of local vertex operators. In other words, the 
partitions must not split the operator O. 

Potential short distance divergences occur in the integral ([4 .41 ) when the mutual dis- 
tances between a subset of the dipoles become small. For simplicity let us assume that 
these are the first M dipole 4'{x2a-i, X2a}, I < a < M < N, and let us contract these 
dipoles towards the origin 0. If we rescale the Xi — Xxi, 1 < i < 2M, the MOPE ensures 
us that when A ^ the non-connected e.v. behaves as 

= A-^^^"Cj...^(01 (/>---(/>)o + 0(A-^^'^" + ^) (4.8) 

" ^ N-M 

M 

where C^...^ is the MOPE coefficient for the contraction of the first M cj) operators toward 
the identity operator 1 and is an homogeneous function of {xi, X2, ■ ■ ■ , X2m-1: X2m} with 
degree —Mdv. The Q(^x-^ d.i^ + terms are associated with the insertion of less relevant 
operators. We can use this result, together with ( [4.ti| ) to evaluate the short-distance be- 
havior of the connected e. v. , by first decomposing it into non-connected e.v. of subsets of 
operators, then using the MOPE and finally re-expressing the result in terms of connected 
e.v. The final result is quite simple. To evaluate the small A behavior of {0(p- ■ ■(/>)o°"" we 
have to consider the restricted partitions of the -|- 1 operators {O, ^, ■ • ■ , (p}, such that 
each subset of the partition contains at least one of the M contracted (p operators, then 
use naively the MOPE for each subset, and then take the connected e.v. for the resulting 
product of operators relative to each subset. The final result is of the form 



(0^->^->);3°"" = A"^^^" (partition) Cj 

n 



restricted 

M N-M partitions 




conn 




(4.9) 



other subsets 
not containing O , _ • i -^V^i 



In ( |4.9|) we have separated the first subset of the partition, which contains O, from the 
other subsets, which do not contain O. Now we simply use the fact that in the contribution 
of each restricted partition all the (■ ■ contain the 1 operator, and that there is at 

least one of the subsets of the partition which contains also some other non-contracted 



operator, for instance the observable O. Now, using for instance (|4.6| ) it is trivial to show 
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that any connected e.v. of products of operators which contains the identity operator 1 
vanishes identicaUy, except of course the single e.v. of the identity operator (1^)'^°"" _ 
This ensures us that each restricted partition in the r.h.s. of ( [4.9| ) gives zero and that in 
fact 

(O (/)■■•(/)(/)■■■(/)) = C'(A-^^^^ + ^) (4.10) 

M N-M 

Now, when integrating over the positions of the 2M endpoints of the contracted dipoles 
hut one (which represents the position of the inserted operator), we obtain that the degree 
of convergence of the configuration is in fact 

w^"" = {2M -1)D - Mdu + D = Me > if e > (4.11) 

This ensures that for e > superficial UV divergences are absent from the e.v. of physical 
observables. 

This result is a strong indication that the perturbation theory for e. v. of physical 
observables is UV finite for e > 0, but is not a rigorous proof. We would need in fact 
to perform a much more detailed analysis of the convergence of the manifold integrals 
when the distances between "nested" subsets of points are contracted toward zero in a 
hierarchical way. This "nest analysis" will not be presented here, but follows the analysis 
presented in the next subsections for the U.V. convergence. 

4-2. Cancellation of Infra-red Divergences 

Before discussing the full renormalization of the model at e = 0, we must check that 
perturbation theory for the infinite membrane makes sense when e — > 0_|_, and does not 
suffer from infra-red divergences. As already discussed in subsec. 2.3., we expect strong 
I.R. divergences to occur for non- vanishing e > 0, but we made the conjecture that at 
a fixed A^-th order in perturbation theory, translationally invariant observables are I.R. 
finite as long as e < D/N (and in the domain < D < 2 where perturbation theory is 
a priori well defined). The arguments developed in the previous subsection to deal with 
U.V. divergences for e > can be adapted to deal with I.R. divergences for e < 0, as we 
now show. 

Let us go back to the manifold integral ( [4.4|) which gives the N-th term of the per- 
turbative series of the e.v. of the observable O 

O^^Hz,} = [ud''x,{0{^^y{{z,})UHx2a-UX2a})r'' (4-12) 
J 1=1 a=l 
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with 

p 

i E qir(zi) 

0{4^}{{z^}) = e (4.13) 

and the additional constraint 

E = (4.14) 

i=l 

Potential I.R. divergences occur in the integral ( |4.12| ) when some of the dipole end-points 
positions Xi become large. We can analyse these I.R. singularities following the general 



method of [p4| , |37[| . The most divergent I.R. singularities are given by configurations when 
a subset of the dipoles goes to infinity. For simplicity, let us consider the case when 
these are the last M dipoles (i){x2a-i, X2a}, N - M + 1 < a < N , with 1 < M < A^. We 
shall prove that when we rescale these dipoles positions by 

x^ px^ ; z>2(A^-M) + l (4.15) 

the integrand in ( [1.1 2| ) scales when p ^ oo as 

{0{^,}{{z,}flt (l>{x2a-UX2a} S HpX2a-l. pX2a})r^ = 0{p-^ ''"-^) (4.16) 

This is sufficient to show that the I.R. degree of divergence of this configuration is 

M {2D - diy) - D = Me - D (4.17) 

To show ( [4.16| ), we first use the fact that the integrand is an homogeneous function 
of the Qi's, the Zi^s and the Xj's, with degree —Ndf under the rescaling 

qi A"'^qi , Zi \zi ; Xj Xxj ; (4.18) 

This implies that the I.R. rescaling of the last M dipoles is equivalent to the U.V. rescaling 
of the observable O and of the N — M first dipoles toward the origin 

N-M N 

{0{qA{{Zi}) n (^{X2a-1, X2a} H (/>{pX2a- 1 , PX2a})S°"" 

a=l a=N-M+l 

= p-'''''' {0{p.^^y{{p-h^})^Yl' c/>{p-^X2a-l,p-'x2a} U </>{x2a-l , X2a})r" 

a=l a=N-M+l 

(4.19) 

To study the large p behavior of the r.h.s. of ( [1.1 9| ) amounts to look at the small distance 
behaviour of O times the product of — M dipoles, when the distances are contracted by 
a factor X = p~^ ^ 0, and when the "external momenta" q^ are simultaneously expanded 
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by a factor X~'^ = . This limit is different from the short distance behavior of products 
of operators studied in subsection 3.2, where the O operator was not involved. However, 
the general method used there to establish the MOPE structure can be easily extended to 
our purpose. Indeed, the only difference is that now we have to study the short distance 
behavior of the product 

N-M 

0{^,}{{Z^}) n (t){x2a-UX2a} (4.20) 
a=l 

under the rescaling 

Qi A"''qi , Zi Xzi , Xj ^ Xxj , A (4.21) 

Since O is a product of exponential operators, similar to those appearing in the Fourier 
representation ( |3.12|) of the dipole operators (p, and since the rescaling (^4.21| ) for the 



external momenta is the same than the rescaling done for the internal momenta hj in 
the MOPE derivation, repeating step by step the arguments of subsection 3.2, we obtain 
that the product of operators ( ^.2C1| ) has also a short distance MOPE expansion, which is 
of the form 



0{ciM{^^}) H^2a-l,X2a} = ^ ^ . . . </,({qa. {^J^ {^.}) " ^(0) (4-22) 

A 



where the sum runs over all local operators A of the form 

i(Eqi)r 

A = monomial(Vr, V r, ■ ■ ■)e * (4.23) 

and where the MOPE coefficient CQ^,,,^{{qi}, {-Zi}, {xj}) is a homogeneous function under 
the rescaling ( [4.21| ), with degree 



degree[C^^.,.^] = -{N - M)i^d - dim[A] (4.24) 

As before this result comes from simple power counting, since dim[0] = and dim[(^] = 
—I'd. The Cq^,,,^ coefficients can be calculated explicitly, but we do not need that at this 
stage. The most singular operator A is 

= e^^^^ ; q = Eq. (4.25) 
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Putting together ( [4.19|) and ( [4.22|) , we obtain the I.R. behavior of the non- connected 

e.v. 



{0<P---<P<P---<P)o = p-'""'Cll__.{V^<P---<P)o + Oip-'""'-'') (4.26) 





We must now use ([4. 61) and repeat the arguments of subsection 4.1 to obtain from ( [4.26| ) 
the I.R. behavior of the corresponding connected e.v. We obtain a result similar to ( [4.9|) : 
we have a restricted sum over the partitions of O and the (p operators, with the condition 
that each subset of the partition must contain at least one of the operators which are not 
I.R. expanded, that is O or one of the first N — M (p. Inside each subset of the partition, we 
perform the U.V. contraction of O (if it is in the subset), and of the first (p, and we obtain 
through the MOPE an operator A, and the connected partition function {A(p ■ ■ ■ (/))o°"". If 
the subset contains O, the most singular A is V^, if the subset does not contains O, the 
most singular A is the identity operator 1. We thus obtain an expansion of the form 

(0</>---</>^--^)r" =p-^''' E "'(partition) c5j,...^(y^(/>---</>)r"x 

(4.27) 




restricted 
partitions 



Mvd-D\ 



other subsets 
not containing O 



Now, when integrating over the positions of the M last dipoles, we see that I.R. divergences 
occur only at e = from the 0{p~^^'^) terms in ( [1.27| ). As long as there is at least one 
subset in the partition without O and with one non-contracted (p, the partition contributes 
with a {l(p- ■ ■(^)o°"" = 0. Thus the only dangerous partitions are those such that all the 
M non-contracted (p are in the same subset as O. These dangerous partitions contribute 
by a factor (Vcjc^ ■ • 

It is at this last step that the invariance of the observable O under global translations of 
the membrane plays a crucial role. Indeed, this means that the total charge q = = 0. 

But from ( [4.25| ) we then have 

Vq=o = 1 (4-28) 
and therefore the contribution of the dangerous partitions is proportional to 

(^ci=6t->)r" = (l^j^)r" = (4.29) 

M 
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and vanishes, since M > 0. 

This estabhshes ( ^.161 ), and shows that "superficial" I.R. divergences are absent from 
the manifold integral ( [4.12| ) when e < D/N, as long as translationally invariant observables 
(with zero external charge) are considered. A complete analysis of the I.R. structure of 
the manifold integrals, taking into account nested I.R. configurations, should be required 
to prove rigorously the I.R. finiteness of perturbation theory for e < 0. We shall not give 
this analysis here, but refer to the general results presented in the next subsections. 

4-3. Superficial U.V. Divergences at e = 

We thus have shown that perturbation theory for invariant observables is UV finite for 
e > 0, and I.R. finite for e < 0. We now discuss in more details the structure of the short 
distance singularities at e = 0, which have been characterized in subsect. 3.3. We then de- 
fine the U.V. subtraction operator for a single singular configuration (divergent molecule). 
This is the first, but essential, ingredient in the construction of the full renormalization 
operator and in the proof of the renormalizability for the model. 

Power counting and the MOPE show that at e = 0, there are two kinds of U.V. singular 
configurations: molecules with one atom, with U.V. degree of convergence w = —D, as 
depicted on fig. ^ (a), and molecules with two atoms, with U.V. degree w = 0, as depicted on 
fig. ^ (b). According to the MOPE, the first case is expected to give UV singularities (when 
the distances within the single atom become small) proportional to the insertion of local 
operators, while the second case will give singularities (when the distances within the two 
atoms independently become small) proportional to the insertion of bi-local operators. We 
shall show that by adding adequate counterterms to the manifold integrands, proportional 
to insertions of these local or bi-local operators, we improve the short-distance behavior of 
these integrands, without spoiling the large distance I.R. convergence. 

In the following we shall denote the integrand in the manifold integral ( ^.4| ) and ( |4.12| ) 

by 

{0{^.}i{^^}) n <y^>{x2„-l,X2J))^" = ^^"""(^.) (4.30) 

a=l 

and we consider it as a function of the 2N dipole end-points positions {xj} only. This is 
justified since the P external points Zi and external momenta are kept fixed, and since 
it appears that no additional U.V. singularities occur when some Xj's are close to a Zi. 
The arguments that we shall present below apply equally to the non-connected e.v. 

N 

{0{<i,}{{z^}) n (/'{x2a-l,X2a}))0 = I{Xj) (4.31) 
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which appear in (|2.7|) , as far as U.V. problems are concerned (but which are I.R. divergent). 
We shall rely on the methods and on the notations used in , where the renormalization 
of the model of a random manifold interacting with a single impurity has been considered. 

• One atom divergences: 

First let us consider the singularity associated to a molecule made of one atom. This 
atom, denoted by V, is made of the 2P end-points of the P dipoles that will be contracted 
to a single point. For simplicity we denote 

\V\ = 2P = number of points in V (4.32) 

Let us consider more closely the behavior of the integrand I{xj) when this contraction is 
performed. For that purpose, we chose a point w inside V, that we call the root of the 
atom V (at that stage the choice of w is completely arbitrary) , and we rescale the distances 
inside the atom P by a factor A 

{Xw + X{xi - Xw) if i eV 
(4.33) 
Xi if i^V 

Under this rescaling, the integrand I{xj) has its small A expansion given by the MOPE, 
and it is of the form 

lix.iX)) = 5^ A'^rH/(x,) (4.34) 

The term of order A*^, T^"'^!, is given by the insertions of local operators A — 
monomial(Vr, V^r, ■ ■ ■) with dimension such that 

a — —Pdiy — dim[A] — —Pdiy+p — qu p > g > integers (4.35) 

and for the non-connected integrand I{xj) it is explicitly given by 

T^^hix,) = C^...^i^j-Xw){A{xJcP---cP)o (4.36) 

A: dim[A]= ✓ 

Only the first terms in the expansion ( [4.34| ), such that 

a<a™ax(P) = -D{\V\-l) (4.37) 
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give an U.V. divergence when the \V\ points in V are contracted towards Xw We can sub- 
tract this divergence by subtracting from I{xj) these U.V. divergent terms, thus obtaining 

IM- E ^^"^^i^j) (4-38) 

but we must be careful with I.R. divergences. The terms with a < —D{\V\ — 1) are I.R. 
convergent, since when aU the Xj in V hut the root x-w go to infinity the corresponding 
integral / H dxjT^^h scale as / rfA/A A'^"'^'-^^-. The term with a = -1^(1^1 - 1) is 

I.R. divergent when e = 0, since the corresponding integral scales as / dX/X. We must 
introduce an I.R. regulator, in order to perform this subtraction as small distances (in 
order to subtract the U.V. divergence), but not at large distances (in order to keep the 
I.R. convergence of the original integral). A simple I.R. regulator consists in subtracting 
the T^'^'^'^'^^ I term only if the distances between all the \V\ points inside the atom V are 
smaller than a fixed distance L. We write the characteristic function 

Xp(xj) = W e{L-\x^-Xj\) (4.39) 

and we thus define the subtraction operation as 

Tpl = ^^"^^ + T^^m..)i . (4.40) 

O" <''"max 

With these definitions, it is now clear that the subtracted integrand 

/ - Tpl (4.41) 



behaves, under the rescaling (|4.33|) , as 



(T rnT)(T (\)) - foix^---'^) h>o X^+OO 



(4.42) 



This ensures both U.V. and I.R. convergence with respect to the integration over the Xj's 
in the atom V, the position of the root Xw, that is of the atom V being kept fixed. 

• Two atom divergences: 

We now consider the divergence given by a molecule made of two atoms, such as those 
depicted on fig. |^ (b). Let us denote Ci and C2 these two atoms, which are disconnected 
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subsets of respectively \Ci\ and IC2I points (note that we may have \Ci \ 7^ IC2I). We shaU 
caU the coUection of these two disconnected atoms 

V = {€1X2} (4.43) 

a diagram, and we caU Ci and C2 its connected components. We now consider the behaviour 
of the integrand I{xj) when both atoms are contracted simultaneously. For that purpose, 
we chose a root x^i inside Ci and a root Xw2 inside C2, and rescale the distances inside Ci 
and C2 by the factor A 

x^^ + X{xi - x^J if i e Ci 
Xi Xi{X) = { x^^ + X{xi - x^^) if z G C2 (4.44) 
Xi otherwise 

Under this rescaling, the integrand I{xj) has a small A expansion of the form ( [4.34|) , 
which is still given by the MOPE, but now in terms of the insertion of bi-local operators 
^{xwn x^u^) with the powers of A given by cr = — (|Ci| + \C2\)diy/2 — dim[$]. Again we 
denote by T'^'^^I the coefficient of the power A*^ in ( |4.34D , and we subtract all the U.V. 
divergent terms such that 

a < a„,ax = -D{\Ci\ - 1 + IC2I - 1) (4.45) 

In fact, power counting shows that for these diagrams with two atoms, there is only one 
divergent term with a = amax, but we may discuss as well the general case where there 
are several divergent terms. Again, the subtraction of the term with a = amax gives 
I.R. divergencies, and we must modify this last subtraction through a I.R. regulator. We 
write the characteristic function for the diagram V as the the product of the characteristic 
functions of its components 

Xt,{xj) = Xc,{xj)Xc,ixj) (4.46) 

and define the subtraction operation 7p for the diagram V corresponding to the divergent 
molecule made of two atoms Ci and C2 by the same formula ( [4.40|) as for one atom. With 
these definitions, the subtracted integrand behaves under the rescaling as 

I-Tpl = o(A-(|Cil+|c^l-2)^) (4.47) 
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fig. 5 

This ensures both U.V. and I.R. convergence with respect to the integration over the Xj 's 
inside the two components Ci and C2 of V. 

• Sum over roots 

We have seen that we need to choose a root, i.e. a contraction point, inside the 
connected components of the divergent diagrams. However we expect that the final result 
for the subtraction operation is independent of the choice of root. For this reason, and 
in order to have a symmetric formulation, it is better to sum over possible roots, with 
an adequate coefficient. For the simple situation discussed here (subtraction of superficial 
divergences), the simplest choice is to choose the same weight for all possible roots, as 
done in [|2^. For the diagram of fig. ^ (a) with one atom V and one root t« G P, we call 
the pair {V^w} a rooted diagram 

= {V,w] (4.48) 

and we assign to each rooting a weight 

W{V^) = ^ (4.49) 

Similarly, the diagram with two atoms Ci and C2 and two roots wi and W2 is called a 
rooted diagram with two components 

= {{Ci,wi},{C2,W2}] = {Ci^^,C2^J (4.50) 

with weight 




(4.51) 



These rooted diagrams are depicted on fig. ^ 

With these notations, since we have by definition 

W{V^) = 1 (4.52) 

roots w 
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fie. 6 



we can replace in ( [4.42| ) and ( |4.47[ ) the subtraction operator Tp, which depends in fact 



exphcitly on the rooting w, by 

Tv ^ W{V^)Tp^ (4.53) 



roots w 



4.4- U.V. divergences and subdivergences: the example of the 2nd order. 

We can construct now a general renormalization operation which makes perturbation 
theory finite. As in standard local field theories, the renormalization procedure is defined 
by an iterative procedure. 

First we subtract the divergent diagrams which appear at the first order in perturba- 
tion theory. These are the diagrams which form 1-loop connected molecules, as depicted on 
fig. Ij. We subtract these 1-loop superficial divergences by applying the subtraction opera- 
tion 7p to all families JF^^) of disjoint 1-loop divergent diagrams JF^^) = {Vi, . . . .Vi}. We 
use the notations of and call two diagrams disjoint if all their connected components 



(atoms) are mutually disjoint. 

We can perform these subtractions by applying to the integrand / the subtraction 
operation 

^^'^ = 1 + E n (-^^)^ (4.54) 

We note that with our definition ( [1.40| ) of the operators 7p, they commute for disjoint 
diagrams Vi and V2- This subtraction is sufficient to make perturbation theory U.V. 
finite at first order (we do not give the proof at this stage). It is simple to see that this 
subtraction operation R*^^) amounts to add counterterms STi^^^ to the bare Hamiltonian 
7i, proportional to the insertion of the operators 1, (Vr)-^ and (j). 

We now want to subtract superficial UV divergences which appear at second order, 
once 1-loop divergences have been subtracted, i.e. for the theory defined by the Hamilto- 
nian Ti + 6H^^\ We thus have to subtract 2-loops divergent diagrams, which are depicted 
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fig. 7 

on fig. 1^ (some diagrams, similar to the third diagram of fig. are in fact not divergent 
and are omitted for simphcity). These subtractions are performed by including in ( |4.54| ) 
families of disjoint l-or-2-loop divergent diagrams. 

But we have also to subtract the new 1-loop divergent diagrams which appear because 
of the additional term 5H^^\ These diagrams are depicted on fig. ||, and contain insertions 
of the operators 1, (Vr)^ and 4), times the corresponding 1-loop counterterms, which 
are associated with 1-loop divergent diagrams. "Expanding" these contractions into the 
corresponding diagrams, we can rewrite the additional subtractions in terms of 2-loop 
divergent diagrams containing 1-loop divergent sub-diagrams, as depicted on fig. |[ 

The corresponding subtractions can be written as products of T-p operators, as in 
( [4.54| ), for l-and-2-loop divergent diagrams, but now some 1-loop diagrams may be "con- 
tained" into 2-loop diagrams (the precise concept of inclusion for diagrams is given in 
and will be recalled below). One thus starts to see the appearance of "nested" subtractions. 

This subtraction procedure can be iterated to 3-loops, etc. . . It produces a renormal- 
ization operation, that is expressed as a sum of subtraction operations over all possible 
families of "nested" divergent diagrams, that we shall call "forests", in analogy with the 
well known Zimmermann's forests of renormalization theory. We now give the general 
formula for the renormalization operation. 

4-5. Renormalization in terms of forests of divergent diagrams: 

In the above sketchy discussion, we also have not taken into account the fact that the 
subtraction operators T are properly defined by ( [1.34|) for rooted diagrams V^, and this 
introduces some subtleties when we iterate the subtractions. To deal with this problem. 



let us recall some definitions of |j2^ 



A diagram "P is a family {Ci,C2, . . .} of (possibly several) disjoint subsets Ci of end- 
points, called the connected components of the diagram. 
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• A diagram Vi is included in a diagram V2 if every connected component of Vi is 
included in a connected component of V2 (not necessarily the same). This is denoted 
Vi -< V2. 

• Let us consider two rooted diagrams Viyj^ and V2W2 such that Vi ^ V2- The rootings 
wi and W2 are said to be compatible if the root of every connected component of V2 
is either a root of a connected component of Vi, or not in any connected component 
of Vi. This is denoted Pi^^^ -< V2w2- 

• The intersection Vi A V2 of two diagrams Vi and V2 is the maximal diagram included 
in Pi and V2 (i.e. such that P ^ Pi and P ^ P2 ^ P ^ Pi A ^2)- 

• Similarly, the union Pi V P2 of two diagrams Pi and P2 is the minimal diagram 
containing Pi and P2 (i.e. such that Pi ^ P and P2 ^ P =^ Pi V P2 ^ P). 
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• A forest JF is a collection of diagrams such that every pair (Vi, V2) of elements of JF 
are either disjoint {Vi A P2 = 0) or included {Vi ^ P2 or V2 ^Vi). The empty set 
is considered as a forest. 

• A compatibly rooted forest J-'q is a collection of rooted diagrams such that the corre- 
sponding diagrams form a forest and the rootings are all mutually compatible (if the 
diagrams are disjoint the rootings are by definition compatible). 

• To any compatibly rooted forest JF^ we associate a weight factor W{J-'^), defined as 
the product over all roots Wi of the inverse of the number of points in the largest 
connected component Ci of the diagrams of JF which has Wi as its root. 

^(•^®) = Ht^ (4-55) 
It is such that one has in particular 

W{:F^) = 1 (4.56) 

compatible 
rootings © 

With these notations, we shall define a renormalization operator R which generalizes 
( [4.54|) and is expressed simply as a sum over all (compatibly rooted) forests of divergent 
diagrams, i.e. of (rooted) diagrams which, together with the dipole which join the end 
points in the various connected components, form a divergent molecule such as those 
depicted on fig. ||. 

The last subtle point in this construction is to define the product of subtractions 
operators Tp-^ and Tp^ for two diagrams such that Vi -< V2- We omit for simplicity the 
rootings Wi and W2, but we have to remember that we assume that these rooting are taken 
into account, and that they are compatible for nested diagrams. The subtraction operator 
for a single diagram V is defined by ( [4.39| ) and ( [4.40| ), i.e. by 



Tpl = y T^'^h + T^'"-^''^/ ■ Xp (4.57) 



O" <''"max 



with I the term of degree cr of the expansion of / with respect to the (rooted) diagram 
V. The last term singles out the "marginal" term, with degree 



components 

= — £>(# of points in P — # of components of V) 
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and X-p is the characteristic function 



= n (4.59) 



components 



which enforce the constraint that the term with degree a"max is subtracted only if all the 
distances within each connected component of V are smaller than the I.R. cut-off L (as 
introduced in (|4.39| )). If we now consider two distinct diagrams V and V in a forest, there 



are two cases: 

- V and V are disjoint {V AV = . In that case, the operators 7p and 7p/ , as defined 
by (|4.57| ), commutes and their action on / is given by 



+ T(^--)T'^'^'V-Xp + T^'^^T'^^'-^^-h -Xj,, (4.60) 

_^ 7-('T„,ax)7-'('^max) J . XpXp> 

V and V are nested. Let us assume that V is included in V' {V -< V). Then two 
operators T^'^^ and T'^'^ which pick the terms of degree a when expanding V and 
the term with degree a' when expanding V' commute too (this is a consequence of 
the MOPE). The only difficult point is to decide how the dilation of V will act on the 
I.R. cut-off Xpi . It turns out that there are different possible choices, which result in 
different renormalization procedures, which are equally consistent. We shall choose 
the most simple: we state that the operators T*^'^-' (which extract the homogeneous 
parts of the integrand) act on the integrands / but do not act on the I.R. cut-offs X. 
This means specifically that il V -< V , the two operators Tp and T-pi still commute 
and their actions on the integrand / is still given by (|4.60 ). Let us note that if we had 



iterated stricto-sensu the subtraction procedure that we have discussed at 1-loop and 
2-loop orders, we would have obtained a slightly different definition, which is more 
complicated to handle. 

This prescription extends trivially to the product of the Tp operators for an arbitrary 
number of diagrams V within a given forest. With this definition, we can now give the 
general form for the subtraction operator R, which acts on the integrand / (of the form 
( [4.30| ) and ( [4.31| )) of the integrals over internal position space. 
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Definition: 

The subtraction operator R is given by tlie sum over all possible compatibly rooted 
forests made of diagrams which correspond to U. V. divergent connected molecules 



R = n (-^^) (4-61) 



Let us stress that in ( |4.61| ) the empty forest (with no elements) is included, and 



contributes by the identity operator 



WmHi-Tp)) = 1 (4.62) 

■Pe0 



which does not change the integrand /. 

We can now state the main result: 
Theorem 

(1) The renormalized manifold integral 



„ 27V 

/ Hd^'xjKIixj) (4.63) 

i=i 



with the integrand I{xi) defined by (|4.31| ), is U.V. finite for e = 0. 
(2) Similarly, the renormalized manifold integral 



„ 2N 

O^I^J = / lld''x,Iir'^--{x,) , (4.64) 

7 = 1 



with the integrand /'^"""(xi) defined by ([4.30|) , is U.V. and I.R finite for e = 0. 



(3) ( [4.64| ) defines in perturbation theory the renormalized expectation value of the ob- 
servable O (given by (|2.6| )), by 

{0>R = f t%;^0if'{..} (4.65) 



N=0 



6ii is here the (new) renormalized coupling constant. The renormalization proce- 
dure defined by the R subtraction operator amounts to add to the bare Hamiltonian 
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Ti/k-QT counterterms proportional to multi-local operators $, with coefficient defined 
in perturbation theory and containing U.V. divergences (i.e. poles in 1/e) of the form 

AH/keT = ^A<i,(6,,e) U (4.66) 

The forest structure of the R subtraction operator and the MOPE ensures point (3), 
i.e. the fact that renormalization amounts to add multifocal counterterms to the Hamil- 
tonian. This ensures that the model of self-avoiding tethered manifolds is perturbatively 
renormalizable in the usual sense. 

In the rest of this section we give the outline of the proof of points (1) and (2), i.e. of 
the finiteness of the renormalized theory. 

4-6. Proof of U.V. convergence 

Subtraction in terms of nests 

First we shall write the subtraction operation R in terms of all possible diagrams, irre- 
spective to the fact that they correspond or not to connected or multiconnected molecules 
or not, and that these molecules are U.V. divergent or not. For that purpose, we recall 
the concept of nest already introduced in , and familiar in perturbative quantum field 
theory. 

• A nestM is a forest with no disjoint elements. In other words, it is a family of included 
diagrams 

M = {v^^V2^-- - ^ Vm} (4.67) 

Again, the empty set is considered as the nest with no elements. 

• Similarly, a compatibly rooted nest A/® is a compatibly rooted forest with no disjoint 
elements, and we associate to it a weight factor VF(A/'©) defined by ( |4.55| ). 



Finally, there is a (purely technical) subtlety associated to diagrams with a connected 
component made of one single point C = {w}. Indeed, subtracting the non-existing 
divergence associated to a single point amounts to perform no subtraction, i.e. to per- 
form the subtraction associated with the empty set 0, and this must not be counted 
twice, otherwise one gets zero. Therefore, as explained in [^, when expressing renor- 



malization in terms of nests, one must consider only the so-called complete diagrams. 
We denote by Q the set of the 2N end-points in the manifold integral ( [4.12| ). A 
complete diagram of ^ is a diagram V which contains all the points in Q. Its connected 
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components Ci may be trivial, if they contain a single point (C = {i}), or non-trivial, 
if they contains at least two points. In other words, a complete diagram is nothing 
but a partition of Q. For instance, to the empty diagram we associate the complete 
diagram V% and the trivial subtraction operator 70 

^ ^0 = {{i}.^(^G} ; ^0=1 (4.68) 

while the maximal complete diagram is 

Vg = {G} (4.69) 

All the definitions of Subsec. 4.5 apply to complete diagrams, and we define a complete 
nest simply as: 

• A complete nest is a nest of complete diagrams of which contains as its smallest 
element 

A/" = {^0 ^^1 ^---^Vm] (4.70) 

Thus the smallest complete nest A/0 is not the empty set 0, but {P0}. 
With these definitions, and those of Subsec. 4.5, we have 

Lemma 

When acting on the integrand I{xj), the operator R defined by the forest formula 
( [4.61|) may be replaced by the equivalent formula, given now by a sum over all possible 



compatibly rooted complete nests made of all possible complete diagrams 

R = -Y,w{:f^) n (-^^) (4-71) 

Proof 

We first have to show that we can replace in (f4.61D the sum over forests made of 



U.V. divergent connected molecules by the sum over forest made of connected molecules 
(U.V. convergent or divergent). This is easy, since if in a forest there is one diagram V 
corresponding to a U.V. convergent molecule, from power counting the action of the T-p 
operator on / gives zero. 

The second step is to show that we can replace the sum over forests of connected 
molecules by the sum over nests. To show this, let us consider the following simple sit- 
uation: Let Vi and V2 be two disjoint diagrams corresponding to two disjoint connected 
molecules, and let us consider the forest made of Vi and V2 

^ = (4.72) 
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It gives in ( [4 .611) the subtraction operator 



Tr.Tr, (4.73) 

From Vi and V2 we can build the union diagram V = Vi M V2 corresponding to the 
(disconnected) union of the two molecules. We associate to this forest JF the three different 
nests 

{V} , {VuV} , {V2,V} (4.74) 



which give in (|4.71| ) the subtraction operator 

-T-p + Tv.T'p + T'p^r-p (4.75) 

The key point is that these three nest are equivalent to the forest JF. Indeed, when applied 
to the integrand /, we have 

{^-rv+Tv,Tv + Tv,Tv)I = (Tr,Tr,)I (4.76) 

This is in fact quite simple. We can rewrite ( [4.76| ) as 

{l-Tp,)Tp + {l-Tp,)Tp + Tp,Tp, = Tp (4.77) 

and use the fact that all these operators probe the scaling behavior of the integrand under 
the dilations of Vi and V2- T-p^ picks the terms with degrees ai < (Ji^s.^ w.r.t. Vi, T-p^ 
picks the terms with degrees (T2 < cr2inax w.r.t. P2, while Tp picks the terms the terms 
with degree a = 0"i + a2 < cr^^ax = crimax + ^'"2max w.r.t. V. Each of the three operators 
on the l.h.s. of (|4.77|) picks one of the disjoint subsets of terms, which are respectively 



and whose union is 



{cTi > (Jiinax) C^l + (72 < Q-lmax + C^2max} 

[02 > 0-2niax, CTl +<^2 < C^lma^ + 0-2max} (4.78) 
{(Ti < Crimax5C^2 < Cr2max} 

{(Jl + (J2 < CTln^ax + Cr2max} (4-79) 



that is the set of terms picked by Tp. Moreover, the I.R. cut-off associated with the "border 
terms" such that ai = crij^ax or 02 = o"2max coincide in the l.h.s and the r.h.s. of ( [4.76| ), 
since by definition (see ( [4.59|) ) Xp = Xp^Xp^. This proves (|4.76| ). 
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This simple example can be easily generalized. To any compatibly rooted nest TV® 
in ( [4 .711) one can associate a unique forest JF^ in ( [4.61|) , made of the connected molecules 



which make the components of the elements of the nest M. Thus to any JF^ we associate 
a class of nests {A/© and these classes form a partition of the set of all nets. 

Then one can show that the sum of the subtractions within a given class reproduces the 
subtraction of the corresponding forest 

E n (-^^) = n (-^^) (4-80) 

Finally, the weight factors W{Mq) of all the nets within a class are equal to the weight 
factor W{J^Q)) of the corresponding forest. 

The third and final step is to complete all the diagrams in each nest by the trivial 
connected components made of the points which are not in the diagram, thus making all 
the diagrams and the nests complete. One easily checks that this amounts to add the 
factor (—1) in front of the subtraction operation, to cancel the subtraction {—Tq^) = —1 
associated with the trivial complete diagram V^. 

This establishes the lemma. 

Decomposition of the integration domain into Hepp sectors 

The next step is to decompose the domain of integration A (over the positions of 
the internal points Xj, or equivalently over the mutual distances \xj — Xj'\) into disjoint 
"sectors" 7i, which generalize the well known concept of Hepp sectors in perturbative 
quantum field theory. Each sector corresponds to a particular subdomain for the distance 
integration, where the distances may vanish in a particular, well defined order. Thus each 
sector will isolate a family of potential U.V. divergences. 

This is done exactly as in p5| . 

• A saturated nest is a maximal complete nest of Q, i.e. a complete nest S made of 
2N + 1 distinct complete diagrams. 

S = {V^^Vi^--- < V2N-1 -< V2N = Vg} (4.81) 

One shows that one goes from a diagram Vr of S to the next one, Vrj^i by fusing two 
distinct connected components of Vr into a single connected component of Vrj^i. 

• To each saturated nest S is associated a extended Hepp sector Ti^ , defined by the 
inequalities between distances inside each diagram Vr of S. Given a diagram V 
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composed of connected components {Ci,C2, ■ • •}, we define the distance dij between 
two components Ci and Cj of V as the minimal Euchdean distance between a point in 
Ci and a point in Cj . 

dij = mm{\x-y\] x e Ci, y & Cj) (4.82) 

As in we associate to V the minimal distance d^in{V) between the connected 
components of V 

d„in(P) = min dij (4.83) 

If V has only one connected component {V = {C}), by definition dminCP) = oo. Then, 
given a saturated nest S defined by ( [4.81| ), the sector H'^ is defined by the subdomain 
of the integration domain A such that 

dminiVn,) < rfn.in(Pl) < ' ' " < rf„in(Pi?) < ' " • < C^min(Pg) = OO (4.84) 



Following p5[, these sectors are disconnected and form a partition of the integration 
domain A 

n^r\H^'=(lS iiS^S' ; \Jn^ = A (4.85) 

5 

Inside a given sector 7^"^, U.V. divergences may occur only if some distances vanish in a 
hierarchical order compatible with the inequalities ( [4.84|) . 

Overlapping divergences 

Then we study the short distance behavior of the renormalized integrand R/ sepa- 
rately within each sector 'H'^ . At that stage, we face the standard problem of overlapping 
divergences: the subtraction operator R contains subtractions relative to diagrams which 
do not correspond to the subsets of distances which are allowed to vanish in a given sector. 

The solution to this problem turns out to be exactly the same as the one of , itself 



inspired from the construction of in perturbative quantum field theory. For a given 
Hepp sector TY"^, we shall classify the nests J\f of diagrams which appear in Eq. ( 4.71j ) 



into disjoint equivalent classes C. Each equivalent class C is uniquely characterized by a 
"minimal nest" A/'min and a "maximal nest" A/'max, and is made of all the nests A/" such 
that A/'min C A/" C A/'max, that is of the nests Af made of all the diagrams of A/'min and of 
some of the diagrams of A/'max- Then we show that it is enough to apply the restricted 



subtraction operator Rc, obtained by reducing ( [4.71| ) to the nests of a given class C, to 
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make the subtracted integrand R^/ U.V. convergent in the sector Ti^ . Summing over aU 
different classes C we get the subtracted integrand R/ and its U.V. convergence in the 
sector Ti^ is ensured. 

This analysis can be repeated inside each sector Ti^ (the decomposition in classes C 
is different), and thus ensures the U.V. finiteness of the renormalized integral. 
We now recall this general construction. 

Let S = {7^° = V% ^ ■■■ ^ 7^2^ = P^} be a saturated nest, and M = {Vq = 

Vq ^ Vi ^ • ■ ■ ^ Vt} be a compatibly rooted complete nest. We denote {w}t the set of 
roots of Vt- For convenience in notations we add to A/® the maximal element Vt+i = Vg. 
Out of the saturated nest S which defines a sector Ti.'^ and of the nest J\f, which defines a 
subtraction operation, we construct the diagrams 

(4.86) 

n'w^^y^Vj = {n'\{Vj\{w}j))wVj 



One shows that 



V^j = Vj ; V^j -< V^/^ . (4.87) 



As in we shall call the set of Vj a tableau. Hence these diagrams form a nest Afmax, 
which depends on S and on A/", and which contains the original nest A/". Note that a large 
number of the elements Vj of Afmax coincide. Moreover, one shows that there is a unique 
minimal nest A/'min such that the above construction gives the same A/max, and that all the 
nests which give A/max contain A/'min- In other words 



AT ^ A/'max A/'min C AT C A/'max (4. 

and all such nests form an equivalent class C5 (A/'min) 

C5 (A/'min) = {AT: ATmin C AT C ACax} (4.89) 

(We do not discuss here the subtleties related to the fact that A/'min and A/'max are rooted 
nests, like Af, and that the roofings of all the nests in an equivalent class Cs must be 
compatible roofings. These points are discussed in details in [p5| .) 



The interest of this construction is that, within the Hepp sector Ti^ ^ the diagrams of 
the minimal nest A/'min are not U. V. divergent, while the diagrams of the maximal nest 
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A/'niax which are not in A/'min are divergent. If we consider the reduced subtraction operator 
Rc relative to an equivalence class of nests 

He = -Yl ^(■^) n (-^^) (4.90) 

Afec veN 

using ( |4.89| ) we can rewrite it as 

Rc = -w{M^^) n (-'^v) n (1-^^) (4.91) 

This makes clear that applying Rc to the integrand /, we retain the counterterms relative 
to the diagrams in A/'min and we subtract the singular part relative to the diagrams in 

Arnax \ Amin • 

To show that this is enough to make the subtracted integrand U.V. finite in the sector 
Ti^ , let us first give a simple example, inspired from the Subsection 7.3 of . 



A simple example 

We take as a "sector" the domain when the distances in a single diagram TZ become 
small. In other words, we consider the nest S with only one non-trivial element. 

s = {7^° = P0,7^^ =7^} (4.92) 

This is not a saturated nest but this is enough for our example. Similarly, we take for 
subtracting nest N a nest with a single non-trivial rooted diagram i.e. we take 

U = {Vo=Vfi,Vi=V} (4.93) 

As a further simplification, we consider that both TZ and V have only one non-trivial 
connected component (i.e. with more than one point), that we denote also TZ and V, and 
we denote hjw the root of V. 

Generically, the integrand subtracted w.r.t. P, (1 — T-p)! , is not U.V. convergent 
when the distances inside TZ become small. If we apply the above general construction to 
our example, two distinct cases must be considered (see fig. |]): 

If w &TZ then we find the tableau 

n = m Po' = {7^^p} 

po = |p| pi = {'JZUV} ^ ^ ^ 
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In both cases the minimal and maximal nests are 



^fmin = M = {V^ -< P?} ; AT^ax = {^o" ^ ^ 'Pi ^ ^1 } (4-96) 
and the restricted subtraction operator Rc for the equivalence class C = {Af; A/'min C Af G 

A^max} is 

Rc = (l-rpi)(-rp)(l-Tpi) (4.97) 

We want to show that when applying R^ to the integrand /, the result is U.V. con- 
vergent when we contract the distances inside TZ. Let us first consider the unsubtracted 
integrand /, and rescale simultaneously the distances inside Pq , V and Vl by the respective 
scaling factors Ai, A and A2. Since these diagrams are nested, we can apply successively 
the MOPE for each rescaling to show that / has a small A expansion of the form 

I ^ Yl ^r^^'^r I^'"^''"''"^) (4.98) 

(Jl,(T,(T2 

Now we apply the restricted subtraction operator Rc to /. Rc-f has an expansion similar to 
( [4.98|) , but from the definition of the T operators, the cr's are restricted by the inequalities 

cri > an,ax(Po) , < f^max(P) , ^2 > Crrm..{Vl) (4.99) 



With amax(P) defined by (^ ). 
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Let us now consider how Rc/ behaves when the distances inside TZ become smaU. 
For that purpose we now perform a similar rescahng of 7?. by a factor At^. Since from the 
MOPE under a rescahng of P by A, / is decomposed into terms relative to the diagram V 
only (the coefficients of the MOPE) , times terms relative to the original diagram where the 
whole diagram V has been replaced by its root w, it is easy to see that under the rescahng 
of TZ by At^, each term /(""i-'^-'^a) of the expansion ( [4 .981 ) of TZcI is homogeneous and scales 
as 

^ j{a^,a,a2) ^j^J^ ^ + ^2 " (4.100) 

Hence Rc/ has an expansion 

Rc/ ^ J^^n" ^^""^^ (4.101) 

and from the inequalities ( [1.99|) 



'-' max (T^o) -^max(P) + ^max {VD (4.102) 

In our example, and have one component, and V\ has one component ilw eTZ, 
but two components if w ^TZ. This implies that 

amax(Po) = -D{\nnV\-l) , a™ax(Po) = -D{\V\-1) , a„,ax(P!) = -D{\nuV\-l or -2) 

(4.103) 

and, since |7^| = |7^ U P| + |7^ n P| - |P|, that 

an > -D{\n\-1) = a^ax(7^) (4.104) 
This last inequality ensures that the subtracted integrand Rc/ is not U.V. singular w.r.t. 

General Argument 

This can be generalized to get the general proof. First let us consider the slightly 
more general case where we have a single diagram TZ and a nest with a single diagram 
V, as above, but now TZ and P are general complete diagrams with an arbitrary number 
of non-trivial components. Out of the nests S and Af defined by ( f4.92| ) and (|4.93|) we 
construct the 4 complete diagrams 

v'l = v vl=Tzy{^}V ^^-^^^^ 
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The minimal and maximal nests, and the subtraction operation Rc are still defined by 
( [4 .961) and ( [4.97|) . A term of respective degrees aj w.r.t Vj has still degree 



= -or" + (7^ - (7? + a} (4.106) 



w.r.t. TZ (this generalizes (|4.100| ), note that aJ] = 0). Hence in R^/ we get only terms of 
degree 

(Tn > -CTn.ax(Po) + '^max(Po) " f^max(P?) + a„ax(P!) (4.107) 



Since we deal with complete diagrams, their number of points is always 2N, and ( 4.58| ) 
reduces to 

o"max(^) = ^(#com(^) - 2A^) ; #com(^) = number of components of V 

(4.108) 

U.V. convergence is ensured if 

an > a„,ax(7^) (4.109) 



A sufficient condition for ( [1.1 09| ) is (since D > 0) 



#com{Vl) - #com(^?) + #com{Vl) > #com(7^) (4.110) 

(we used #com('Po) = #com(^0) = 2A^) or equivalently 

#com(7^ V{^} V) + #com(7^ AV)> #com(^) + #com(7^) (4.111) 

( [1.111|) follows from the fact that by construction 

7^ V{^} V <nyv #com(7^ V{^} v) > #com(7^ v v) (4.112) 

and from the general inequality valid for any partition A and B of a finite set Q (remember 
that complete diagrams are nothing but partitions of ^) 

#eom {AyB)+ #com {A A B) > #eom (A) + #com (S) (4. 113) 

q.e.d. 

We now consider the general case. Let 

5 = {7^° = P0 ^ ^ — < 7^^^-^ ^ n^^ = Vg} (4.ii4) 
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be a saturated nest and 



Afmin = {Vo = V^<Vi<---^Vt} (4.115) 

be a minimal nest (w.r.t. S). To ATmin is associated the maximal nest 

Ar„,ax = {Pj; / = 0,---,2iV, J = 0,---,T} (4.116) 

and the subtraction operator Rc defined by ( |4.91|) . Using the factorization properties 
which follow from the MOPE, a term of respective degrees Uj w.r.t. Vj has degree 

^7^^ = (4.117) 
J 

w.r.t. IZ^ . U.V. convergence in the Hepp sector Ti^ will be ensured if we have 

a7^I > a^ax(7^^) , / = 0,---,2iV (4.118) 

Since in ( ^4.117| ) the = Vj G A/'min, and the Vj G A/'max, in the expansion of Rc/ w.r.t. 
7?-^, we get only terms with degree 

(^n' > $^^^max(7^i)-CX^ax(^j) (4.119) 
J 

Hence a sufficient condition for ( fl.ll8|) , i.e. for U.V. convergence, is that for all I we have 

$^a™ax(7'j) -amax(^j) > a„ax(7^0 (4.120) 
J 

or equivalently, since we deal with complete diagrams, 

T 

J2 #com(P5) - #com(Pj) > #con.(7^) " 2N (4.121) 
J=0 

We thus have reduced the problem to a given /, i.e. to a given element of the saturated 
nest iS, TZ^ , but we still have an inequality involving all the J's, i.e. all the elements Vj of 
A/lnin- We now show that we can reduce the problem to a simple one inside each Vj. Let 
us call 7^(j_|_j^^ the restriction of 71^ inside "Pj+i 

7^^J+l) = AVj+i (4.122) 
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The diagram Vj, defined by ( [4 .861) , depends on TZ^ only through its restriction 7^(j_|_]^) to 
Vj+i, since one has 

^5 = ^(J+i) V{.},, (4.123) 



And noting that from (|4.122 



7^^o) = P0 , 7^^^+l) = (4.124) 

we may rewrite the r.h.s. of ( |4.121D as 



#eom(7^) -2N=J2 #com(7^fJ+l)) - #com(7^fJ)) (4.125) 
J=0 

Hence, a sufficient condition for U.V. convergence w.r.t. TZ^ is that for each J, one has 

#com(Pj) - #com(Pj) > #com (7^f j+l) ) " #com(7^[J)) (4.126) 

But we can now check this inequahty inside each component of Vj+i separately. Indeed, 
considering each component of Vj+i as a whole Q, the restriction of Vj to this component 
is a complete rooted diagram V of with roots {w}, the restriction of 'R.f^j_^_i^ is another 
complete diagram TZ, the restriction of "^[j) is nothing but TZ AV and the restriction of 
Vj is nothing but 7Z V^w} Thus the inequality ( fl.l26|) reduces to 

#com(7^ V{„} V) - #com(P) > #com(7^) " #com(7^ A V) (4.127) 

This is nothing but ( [4.111|) that we proved above. 

We have thus proved that within a Hepp sector Ti^ , each restricted subtraction oper- 
ator Rc makes the integrand U.V. convergent. Since 

R = ^Rc (4.128) 

c 

the full subtraction operator R makes the integrand U.V. convergent in all Hepp sectors, 
hence in the whole integration domain A. This ends the proof of U.V. convergence, that 
is of point (1) of the theorem of Sec. 4.5. 
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^.7. I.R. finiteness and other points: 

It remains to show point (2) of the Theorem, namely that applying the same sub- 
traction procedure to the connected amplitudes defines renormalized expectation values of 
invariant observables which are both U.V. and I.R. finite at e = 0. 

It is not difficult to prove U.V. finiteness, along the same line as above. The connected 
manifold integrand /'^°""(xj) in ( [4.64|) , defined by the connected e.v. (|4.30|) , has similar 



factorization properties under the short distance expansion as the manifold integrand I{xj). 
Indeed, they can be written as sum of products of /(xj)'s, and applying the MOPE to 
amounts to apply the MOPE to the /'s. Then we can repeat step by step the arguments 
of Sec. 4.6, to prove U.V. finiteness of the renormalized manifold integrals ( [4.64|) . 



Similarly we can show that renormalizing the connected manifold integrals amounts 
to add to the bare Hamiltonian the same counterterms as when renormalizing the non- 
connected manifold integrals. 

The most important point is to show that the connected manifold integrals ( [4.64| ) 
are also I.R. finite, i.e. that no additional divergences occur when some of the distances 
become very large {xj oo). We shall not give a detailed proof of this result here, 
since it is lengthy and not particularly illuminating. We have to follow the procedure of 
Subsec. 4.2, which was used to show I.R. finiteness of unrenormalized connected manifold 
integrals, checking that the U.V. counterterms do not introduce additional I.R. divergent 
terms. This is a consequence of our definition of the subtraction operation Tp, given by 

T-pI = ^ T('^)/ + T^'^"-'')/ ■ X-p (4.129) 

O" <fmax 

When we dilate the distances inside the diagram P by a global factor p, the first terms 
T^^^I scale as p"', with a < amax, and are I.R. convergent when p — > oo (since they are 
U.V. divergent when p ^ 0). The only dangerous term is the last one, with a = cTmax, 
which is logarithmically divergent at e = 0. However, because of the I.R. cut-off xv 
(defined by ( [4.39|) ), this term vanishes if some of the distances inside V becomes large, and 
is not present in the large p asymptotic expansion of T-pI. This is the reason why this I.R. 
cut-off was introduced. This implies that the U.V. counterterms are not I.R. divergent at 
e = 0, and since the bare connected manifold integrals are I.R. finite, that the renormalized 
manifold integrals are also I.R. finite, as long as the external points Zi (in the observables 
0{^.].({zi})) stay at a finite distance from each other. 
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Another subtle point, that we have not discussed, is related to the distributional 
character of the manifold integrals. We have checked that the U.V. power counting of the 
subtracted integrands R/ is good enough. However, if we had dealt with functions, which 
had to be integrated with a well defined measure d^Xj, it would be necessary to check that 
the U.V. bounds are uniform over the integration domain A. This point has been checked 
explicitly for the simpler model of describing a membrane interacting with a single 
point. It would be possible to study this problem for our model, although the techniques 
of ||2^ become very heavy. However, one has to remember that the measure d^Xj is in fact 
defined as a distribution over distance space. This measure is singular on the boundaries 
of A, and these singularities are treated by a finite part prescription (i.e. by integration by 
parts). To have a complete proof of U.V. finiteness, one should check explicitly that these 
integrations by parts do not interfere with the U.V. singularities of the integrand (which 
occur on a much smaller subset of the boundary of A). Explicit calculations at first and 
second orders show that this is indeed the case, and we do not see any reason why this 
might not be the case at higher orders, but we do not have a full proof of this claim. 



5. Scaling Behaviour 
5.1. RG Equation 

The outcome of the previous section is that the theory is renormalizable through ex- 
actly two renormalizations, a wave function renormalization and a coupling constant renor- 
malization. More precisely, let us define a renormalized Hamiltonian Ti.^ [r r] , function of a 
renormalized field rji{x), of a renormalized coupling constant 6r. and of a renormalization 
momentum scale /x, by 

^^[?.]/kBT=^ Jd^x (V.f,(x))V^^^^|^ Jd^x jd^x' 5\v^{x)-v^{x')). 

(5.1) 

One can build in perturbation theory two renormalization factors Zihj^) and ^5(61^), which 
are functions of 6^ and of e, so that the theory, when expressed in terms of the renormalized 
quantities and h^, is UV finite when e 0. This means that if we consider any physical 
observable O(rii), such as the two point function 

F = {v^{x)-v^{y)f , (5.2) 
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its average value 



is a UV finite function of fe^ as e — * 0, when evaluated with the renormalized Hamiltonian 
H^. As usual, the renormalization factors Z and Zb, when expanded in b^, have UV 
divergences, which are poles in e 



Z{b^) = 1 + + + 0{b, 

Zb{b^) = 1 + b^i- + ...) + Oib 



(5.4) 



A major consequence of this renormalizability is the existence of renormalization group 
equations, which as usual come from the freedom in the choice of the renormalization scale 
^, and which can be derived as follows. The renormalized Hamiltonian ?i^[rii] is equal 
to the bare Hamiltonian 7i[r] given by (|2.1| ) by the change from renormalized to bare 
quantities 

r = Z^'\b^)r^ ; b = b^ Zi,{b^) Z''/\b^) fi' (5.5) 
This means for instance that 



{{rM)-rM?)^ = Z-'{{v{x)-v{y)f) . 



(5.6) 



Considering that b is fixed and varying |U, (5^) defines a renormalization group flow for 
6r, characterized by the Wilson [3 function 



d_ 



(5.7) 



Similarly, the anomalous dimension for the field is given by 

d 



p{b^) = li 



ln(rR,/x'') 



1 d 



(5.8) 



InZ 



where u = (2 — D)/2 is the canonical dimension of r. Since the functions P and u are 
dimensionless, they depend only on b^ (and on e and D), and not on fj,. A variation of 
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the renormalization scale n for the two-point function ( ^.31 ) is absorbed into a change of 
6r. and a wave-function renormahzation, through the renormahzation group equation 



F^{x-y) = 



(5.9) 



From (5^) one derives the CaUan-Symanzik equation, which gives the variation of the 
two-point function under global rescaling of the distance £ = |x — y|. Thanks to the simple 
homogeneity relation 

= (5.10) 

the renormalization group equation gives 



^|+/3(6.)A_2.(6.) 



. 



(5.11) 



As usual, if there is an infra-red stable fixed point h^* > 0, such that /9(6r*) = 0, (3'{b^*) > 
0, the large £ behaviour of -Fr(^) will be governed by this fixed point, leading to 



V 



(5.12) 



(5.13) 



i.e. to the fractal dimension dfrac of the membrane 
5.2. 1-Loop Results 

We now apply the MOPE formalism to compute the renormalization group functions 
for the model of self-avoiding tethered membranes at one loop, and show how to recover 
the results of [1T2|-[T3-18| for the scaling indices at first order in e. For simplicity let us 
denote by (f>{x, y} the two-point operator 



(pix^y} = S'^irix) -r{y)) 



(5.14) 



The U.V. divergent diagrams at first order are represented on fig. |10 




X2\^ — W3^2 



X3 



Xj X2 



yi 



xi 




X2 



fig. 10 
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We need the MOPE for the contraction (/}{xi,X2} local operators when xi X2, 
which is 



4>{xi,X2} = Ci{xi - X2} + C"^{xi - X2} :VQr(x)V/3r(a;): + . . . (5.15) 
with X = {xi + X2)/2. We start from the OPE 

1 + i(x^ -x^):Var'^(x):k^ 



gik(r(xi)-r(x2)) 



gk2G(xi-X2) 



(5.16) 



The MOPE ( ^.15|) is obtained by integrating over k , leading to 

Cj{xi-X2} = (47r)-'^/2(-G(xi-X2))-"/' (5.17) 

C'/ixi - X2} = - ^ (47r)-"/2 (-G(xi - X2))-^-"/'(x? - x^)(x? - x?) (5.18) 
The MOPE for (jxp ^ (p is of the form 

H^i, yi}<P{x2, y2}^=^^ cj^{xi - X2, 2/1 - 2/2} H^, y} + ■■■ (5.19) 

HI — >H2 

where x = (xi + X2)/2, y = {yi + y2)/2. We use the OPE 

giki(r(xi)-r(yi))gik2(r(a;2)-r(y2)) _ gi(kir(xi)+k2r(x2))g-i(kir(yi)+k2r(y2)) 

^ J^gi(ki+k2)F(x) _^ j g-kik2G(xi-x2) j^g-i(ki+k2)?(y) _^ j g-kik2G(yi-y2) 

(5.20) 

The coefficient C^^ is obtained by integrating over ki and £2, after inserting the identity 
1 = / d'^k 5'^ (k. — ki — k2) so that we get for (pep 

I If ^ ^ d^k5^(k - ki - k2) e''^(''(^)-''(^)) ^-^.^.{G{.,-.,)+G{y,-y,)) (g^l) 

J J J {2txY (27r)^ 

After expanding 5^(k — ki — k2) = 5^(ki +^2) + ■ ■ ■, the integral over k builds the operator 
(/){x,j/}, and the integral over ki gives the C^^ coefficient 

C^^xi - X2, 2/1 - 2/2} = (47r)-'^/2 [-G(xi - X2) - G{y^ - 2/2)]"'/' (5.22) 
Finally, the MOPE for the one-particle-reducible diagram is easily found to factorize 



into 



(/){xi,yi}^{x2,X3} = Ci{x2,X3}(^{xi,yi} ■ • • (5.23) 
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This shows that the divergence naively associated to this diagram is in fact not present, 
and is a reflect of the relevant divergence (p ^ 1. 

Let us now compute the renormalization factors Z and Zf, at flrst order in b^, i.e. the 
coefficients A and B in ( |5.4|) . One starts from the action ( |5.1| ) without counterterms, that 
is with Z = Z\j = \. The perturbative expansion gives at order h^^ a term 

/ / ^^''^it^''^2</'{xi,X2} (5.24) 

When x\—^X2, it gives a singular contribution (a pole in e) equal to 

[d^x :(Vr,(x))2: / - X2) ^ Cf{x^ - x^] (5.25) 

^ J J\xi-X2\<L 

L is an arbitrary I.R. length scale, and since the residue of the pole does not depend on L, 
we can set it to L = (the inverse of the renormalization mass scale). The last integral 
over (xi — X2) can be performed explicitly, and gives 

[ d^'y - Cr{y} = ^ — (47r)-^ (2 - D)^+^ ( ^^^] ' (5.26) 
This pole is therefore subtracted by a wave-function renormalization, by choosing 



At order we get a factor 



l + ^ji^(4.r«(2-I,r»^^j (5.27) 



''hV' " " d0^.jdOj^jdDj,^<i%^{i.j,„j}^{i.2,j,2} (5.28) 




When xi — > X2, yi ^2/2 simultaneously, it gives a singular contribution equal to 

Jj d^xd^y<P{x,y} y"^^^_^^d^(xi-X2)d^(yi-y2)Cj4xi-X2,2/i-2/2} (5.29) 
The singular part of this last double integral can be calculated explicitly, leading to 



1 (2-^r^ / 27r^ ^(2^ 
e (47r)l Vr(I)/2)J p 



+ 

2^ 



(5.30) 
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The same pole comes also from the situation where xi ^ y2, yi X2- These two poles 
are subtracted by a coupling constant renormalization, by choosing 

+ f ( 



2 r 



D 
2-D 



2D 
2-D 



(5.31) 



Prom ( |5.27| ) and (|5.31| ) is is straightforward to compute the renormalization group 
functions, using (p3|), (O) and (|5.8|). We find 



bj (2-D)-^+i / Y 
2 (47r)l Vr(I)/2) J 



rf(2-D)2 
P , ^ 2 



. 2 - D 6, {2-D 



2-D 



2n^ \ 



2+i 



(5.32) 



(5.33) 



\T{D/2)J 

The coefficient of P is positive for < -D < 2. Therefore for small positive e > 0, 
the /3 function has an I.R. attractive fixed point b^* oc e, which governs the large distance 



scaling behavior of the membrane. Using (|5.12| ) the anomalous dimension at order e is 



D 



V 



+ e 



4D 



D 
2-D 



(2-Dy r ( 

^ \2-D 



(5.34) 



This result agrees with that given in |]T^JT^,|TB[ (it differs only by a term of order 0{e 



and one has to remember that the expansion parameter e = 4D — d{2 — D) used in these 
papers is not the expansion parameter e used here, but is £ = 2e). 



6. Finite Membranes 

6.1. Extension of the model to finite membranes 



Up to now we have considered infinite membranes only: the model defined by the 
action (|2.1| ) is defined on the infinite Euclidean space IR^. It is important to extend 
our results to finite membranes. In particular, the direct renormalization method used 



m 



r2| , p!3|JT^] to perform one-loop calculations makes a crucial use of the finiteness of the 
membrane to define renormalized quantities, such as the dimensionless second virial co- 
efficient g, and to study its fiow with the size of the membrane z as z ^ oo. In order 
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to prove the consistency of the direct renormahzation method, we need to relate it to the 
field theoretical renormahzation method developed here. Moreover, several scaling indices 
characterize finite membranes. This is the case of the configuration exponent 7, and of the 
contact exponents ^'s for membranes with boundaries. 

Let us first consider a finite tethered membrane without boundary. It corresponds for 
instance to an irregular lattice with coordination defects and link varying elastic moduli. 
Since on such a lattice there is no natural coordinate system, it is better to describe such 
a system in the continuum limit in a coordinate independent way. This means that we 
choose arbitrarily some coordinate system {x} on our manifold Ai (this is always possible 
locally, and we want to find a continuum Hamiltonian which describes the elastic properties 
of such a membrane. If we first neglect self-avoidance and consider a phantom membrane, 
its elastic properties are simply encoded by a fixed internal Riemannian metric gctp{x) on 
the manifold Ai. The elastic Hamiltonian is still Gaussian, and of the form 

^ei[i^/kBT =11 d''x^^d^v{x)g"f'{x)dpv{x) (6.1) 
^ Jm 

As before, a = 1, ■ ■ ■ , -D label the internal coordinate indices. This Hamiltonian transforms 
covariantly under changes of coordinates: the same membrane, with the same internal 
elastic properties, is described in the coordinates x' by the metric g'^p{x') = ■§^-§^g-ysix). 

It is now easy to generalize the Hamiltonian (|2.1| ) for self-avoiding membranes. The 
covariance under a change of coordinates fixes it to be 

n[r\/]^^T =11 d''x.Mx)D^v{x)Dj{x) 

r r ^'-'^ 

+ - d'^x^f^ d'^x' ^^5\v{x)-v{x')) 
^ JM Jm 

The Da and are covariant and contravariant derivatives for the internal metric g. 
Formally, the perturbative expansion can be constructed along the same lines as for an 
infinite membrane. It involves the same diagrams and the same manifold integrals ( 2.15| ) 
and ( p.lTp . Simply, the two point function G{x, y) on IR^ given by (|2.13|) has to be replaced 
by the massless propagator on the surface ^A, GM{x,y), solution of 

-A.G«(x,,) = -J=*>-,)-^ (6.3) 

where A = D^Da is the scalar Laplacian on A4 with respect to x. The additional source 
term, proportional to the inverse of the area of the manifold Vol(A^) = d^x \Jg{x), 
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is here to take care of the zero mode of the Laplacian. This zero mode is a constant 
and is associated with the global translation invariance r(x) r(x) + Fq of the model. 
The propagator, solution of ( |6.3|) , is defined up to a constant. However, this arbitrariness 
disappear in the partition function ( p.3|) and in the translationally invariant observables 
( p.7|) . Indeed, the determinant A{xi} and the minors N^^{xi^ Zn} involve only differences 
of propagators Gm{x, y). 

For open membranes, i.e. membranes with a boundary, one must specify the bound- 
ary conditions. Translation invariance implies that the natural boundary conditions are 
Neumann boundary conditions: the normal derivative D±r of the field r must vanish on 
the [D — 1) -dimensional boundary dAi of A4 

D±r{x) = n''{x)Dar{x) =0 if x E dM , (6.4) 

where n" is the unit vector normal to the boundary. 



6.2. Examples of D- dimensional manifolds 

The above described construction is somewhat formal. In practice perturbative cal- 
culations have to be performed for e = 2D — {2 — D)d/2 small and non-integer < -D < 2. 
A general formulation of curved (Riemannian) manifolds with non-integer dimensions is 
lacking, but one can device several explicit examples where calculations are possible. 

A simple case is discussed in |^5|: the constant curvature (i-dimensional sphere So 
with radius L. By embedding the sphere in IR"^^^, one can generalize the concept of 
distance geometry in IR"^ to So, and use 80(1) -|- 1) invariance (instead of Euclidean in- 

M 

variance E(L>)) to write the measure n d Xi^J g{xi) for M points over 5^) as a distribution 

i=l 

d^Mi^iji L} over the M{M — l)/2 geodesic distances iij between these points. This 
distribution has a simple form when expressed in terms of the so-called chord distance on 
the sphere (see for details). This distribution is a meromorphic function of D, and 
allows for an analytic continuation to < D < 2. The massless propagator ^^^(x, y) can 
also be defined on Sd for non-integer D: it is a function of the geodesic distance i between 
X and y only (thanks to SO{D + 1) invariance), and has for instance the following integral 
representation 



GsM = L'-""^^^^ T- hl + r2-2rcos(£/L)|^ -|l-r|^-^l , (6.5) 



47r 2 Jo ^ 
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where L is the radius of the sphere, which aUows for an exphcit analytic continuation to 
< D < 2. Eq. ( |6.5| ) can be derived simply by starting from the massless propagator in 
infinite flat space in D + 1 dimensions. Details are given in appendix B. 

A more general class of examples is provided by taking for the manifold Ai the product 
of a "physical" 2-dimensional manifold and of a (-D — 2) -dimensional sphere Sd-2- 

M = M2®Sd-2 (6.6) 

(2) 

M.2 has a fixed metric (7^^, and Sd-2 has the constant curvature metric of the sphere with 
radius L. When L oo we get M = M2® 1R^"^ Points X on M. are labeled by their 
coordinates x in M.2 and x' on Sd-2, x = {x,x'). The integration measure for P points 
on M is the product of the ordinary measure Ild'^x^/g^'^\x) on M.2 and of the distance 
measure dfj,p{i'^j) on Sd-2- The massless propagator Gjiiix^y) on M is easily obtained 
from the heat-kernel on A^2, -ft^A^2(^5 2/5 ^) that on Sd-2, Ks^^^^i^''^ where I' is the 
geodesic distance between x' and y' , by the relation 

Gjiiix^y) = / dsKM2ix,y,s)Ksi,_2{i';s) . (6.7) 
^0 

Similar tricks can be used to define the model on open manifolds, for instance the D- 
dimensional ball Bd- 

6.3. MOPE and UV Singularities for finite Membranes 

Within this framework, the analysis of the short distance singularities of the pertur- 
bative expansion of the model can be performed along similar lines as for the infinite fiat 
manifold. In that case, the MOPE (|3.25| ) was derived from the OPE for the product of 
exponentials of the r field inside the atoms P, given by ( p.l3|) , (p. 141) and ( |3.15| ). A similar 



OPE holds in curved space |^^: when the points Xi inside an atom V tend toward the 
point x-p, one has 

- r ^ -| — 5 S k.i-iijGo{yi—yj) 

Here yi — Xi — x-p and k-p = "Eiev ki is the total momentum of the atom V. Go{yi, yj) is 
the free massless propagator in the infinite fiat space with the metric gctp{xp) at the point 
xp: 

Go{y^,yj) = — ^{2-d)Sd ' \\yi~y^\'^^v = (y? -yj)iyf -yj)9»f3ixv) (6.9) 
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There are two crucial differences between ( |6.8|) and the OPE in flat infinite space . 

(1) The sum over the operators Ai^x-p) is now performed on operators which are not only 
monomials of the derivatives of r (Vr(x-p), V^r(x-p), . . .), but depend also on the 
metric tensor g and of its derivatives {g{x'p), Vg{x-p), V'^g{x'p), . . .) at the point x-p 
(V = d / dx is the ordinary derivative). 

(2) The standard normal product in flat space :[ ]: (deflned by factorizing out the ex- 
pectation value of the single operator [ ], i.e. by subtracting all tadpole diagrams 
in the perturbative expansion of the expectation values of observables involving the 
operator [ ]) is replaced by the general normal ordered prescription in curved space, 
A/'[ ], which is deflned in Appendix C. The basic idea is to make the operators U.V. 
flnite (in curved space) by subtractions which do not depend explicitly on the back- 
ground curvature. This prescription coincides in flat inflnite space with the usual : [ ] : 
prescription. 

The coefficients of the OPE are still monomials in the and are homogeneous 
functions of the y^'s. They depend only on the local metric g at the origin x-p, since 
they are polynomials of the \\yi — and in particular they do not depend on the 

derivatives of the metric at the point x-p. A direct consequence is that for the operators A 
which do not involve derivatives of the metric, that is for the operators which are already 
present in the OPE in flat Euclidean space, the coefficients C"^{yi,ki} are the same as 
those in ( |3T5| ). 

The outline of the derivation of (|6.8| ) is given in Appendix C. It relies on the short 
distance behavior of the propagator G{x, y) in general curved space, which itself can be 
written as an OPE involving local operators A which depend on the metric g and its 
derivatives. 

Following the derivation of the MOPE ( |3.25| ) given in Section 3.2, we can now derive 
a general MOPE in curved space. The multilocal operators are of the form (similar to 
(PI)) 



{A-p,m-p 



y{xp} 



d'^r U <|A/' 
veM 



Arixp){iVrr^ 5"(r-r(xp)) 



(6.10) 



where the A-p are the general local operators depending not only on derivatives of r but 
also on the metric and its derivatives. Let us consider a product of M such multilocal 
operators let us partition the set of points {x} into a set of atoms {V} and let the 
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points X which belong to each atom V collapse toward x-p. The short distance behavior of 
the product is then given by a MOPE 

n = Y.^t.:'^M{y}-H^v} (6.11) 

where as usual the y's are the relative distances inside each atom V {y = x — x-p). The 
coefficients C^-^^ ,,,^^^^{y} depend on the metric g only through the \\yi — j/jH^^'s inside 
each atom V, and therefore on the metric g{x-p) but not on the derivatives of the metric 
at the xp's. These coefficients are homogeneous functions of the y's, whose degree can 
be extracted by power-counting analysis from the dimensions of the operators A and of 
the m's in the multilocal operators As for the OPE, is is important to notice that the 
coefficients of the MOPE C*^ ... ^^^^{y} for operators $'s which do not involve derivatives 
of the metric are the same as those of the MOPE in flat Euclidean space. 

As a consequence, the short-distance UV divergences will be proportional to the in- 
sertion of relevant multilocal operators of the form ( |6.10| ), i.e. operators such that /■ ■ ■/ $ 
has a scaling dimension > at the critical dimension (corresponding to e = 0). Moreover, 
since the dimensional regularization prescription that we use is known to preserve the 
covariance under changes of coordinates x x'{x) on the manifold the U.V. diver- 
gences (i.e. poles at e = 0) will be proportional to multilocal operators which transform 
covariantly under such changes of coordinates. 

6.4- Renormalization for finite Membranes 

Following the arguments of section 4, the theory on a curved space will be renormalized 
and made U.V. finite by adding the adequate counterterms to the Hamiltonian (|6.2| ), which 
are proportional to all relevant and marginal multilocal operators. Thus, the situation for 
the renormalization of the model of self-avoiding membranes with a fixed internal curved 
metric is similar to that of local field theories in curved space: new counterterms have to 
be added, corresponding to local, metric dependent, operators. 

For a closed manifold A4 and for dimension D < 2, the only marginally relevant 
operators are those which are in the Hamiltonian (^]5). The renormalized Hamiltonian 
which makes the theory finite for e = is therefore 

-Hn/kBT =^M. f d^xy^D''r^{x)D^r^{x) 

+ lb^lJ.'Ztib^) [ d^xv^/ d^x'vWy5"(?K(x)-r^(x')) , 
^ Jm Jm 

(6.12) 
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where ^(6r) and are the same renormahzation factors as those calculated for the 

infinite flat space model (see ([O]), ( J5.4D ). Operators which depend on the derivatives of 
the metric, that is in a covariant form on the curvature tensor and its covariant derivatives, 
are irrelevant at e = for D < 2. The less irrelevant one is just the integral of the scalar 
curvature R{x), 



d'^xy/g{x)R{x) , (6.13) 

M 

which has scaling dimension 2 — D. Other operators, such as J ^Jg R°'^DctrDpr^ 
f ^/g R f y/gd'^{r — r), etc. . . have larger scaling dimensions. 

For open membranes with a boundary, additional boundary operators have to be taken 
into account. Indeed, as for local field theories, the short distance behaviour of products 
of multilocal operators is modified when some subset of points collapse toward a point 
on the boundary dAi of the manifold A4. The MOPE (|6.11| ) now contains multilocal 



operators of the form (|6.10|) with special boundary local operators B{z), z G in 



addition to the bulk operators A{x). Additional short distance divergences may occur on 
the boundary and have to be subtracted by additional counterterms proportional tothe 
relevant multilocal operators containing boundary local operators. 

For D < 2, the most relevant boundary operator is simply the boundary integral 



d^'-'y^Mj), (6.14) 

dM 

{h[y) being the metric on the boundary induced by the metric g in M). At first order, it 
is generated when the two end-points of a dipole collapse toward the boundary. Its scaling 
dimension is equal to 1 — D, and since L is a purely geometrical term, its expectation 
value and its scaling dimension are not modified by the fluctuations of r. Thus it is 
irrelevant for D > 1, marginally relevant for D = 1, and strongly relevant for D > 
1, like the volume integral jj^d^x \/ g{x) for D > 0. A higher dimensional operator 
is Jgj^ \/h{y) k{y), with k{y) the extrinsic curvature of the boundary in Ai, and has 
dimension 2 — D, it is irrelevant for D < 2. Other irrelevant operators are for instance 

IdM VW) Im - ?(^)), IdM VW) loM VWlS^im - ?(?/')), etc. . . 

This implies that for D = 1, i.e. for the case of polymers, an additional counterterm 
of the form 

Zdib^) [ d^'-^^Mj) (6.15) 

JdM 
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has to be introduced to make the renormahzed theory finite. The renormahzation factor 
Zq depends only on 6pt7 B-nd. CB/D. be viewed. 3,s 3; "boundary coupUng-constant" which has 
to be added to the Hamiltonian to cancel the (divergent) boundary term generated by 
the fiuctuations of the membrane. This new term does not modify the already existing 
renormalizations of the bulk operators, Z{h-^ and Z{,(6pi). This renormahzation factor is of 
course not new, and can be recovered from the standard treatment of the Edwards model. 
If L> 7^ 1, the boundary operator is not superficially relevant, and no renormahzation is 
required if we use dimensional regularization. The renormahzed Hamiltonian ( |6.12|) is 
sufficient to make the theory UV finite for e = 0, both for closed and for open membranes. 

6.5. Finite Size Scaling and Hyperscaling 

The fact that the renormahzation factors are the same for the theory in infinite D- 
dimensional space and in a finite space is very similar to what occurs for local field theories. 
This property of local field theories allows to derive the so-called finite size scaling laws 
for critical systems in finite geometries. In the same way this property of self-avoiding 
membranes allows to derive scaling laws for the behaviour of large finite (closed or open) 
membranes. 

As a first application, let us study the scaling behavior of the partition function of a 
closed membrane Ai. The partition function is defined in the standard way as the sum 
over all configurations r, with only one difference with respect to (|2.2|) : in order to remove 
the translational zero mode, which gives an infinite factor equal to the volume of the 
(i-dimensional space in which the membrane fiuctuates, one must fix the position of the 
center-of-mass rem 



with the Hamiltonian 7i[r]/kBT given by ( |6.2| ). 

If we now scale the internal size of the manifold 7W by a factor X , for instance by 
rescaling the metric ga.fi{x) X'^gap^x), we expect that in the large X limit, the partition 
function Z(h,X) of the new manifold scales as 





(6.16) 



The partition function is then 




(6.17) 



Z{h-X) = X^-i (cst.)^°^(-^)- 



(6.18) 
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where est. is a non universal constant (depending for instance on the way the model is 
regularized), but independent of the shape and the size of the manifold, and 7 a univer- 
sal (independent of the regularization) configuration exponent. As we show now, if the 
dimension D of the membrane is not an integer (L> 7^ 1, 2) this configuration exponent is 
simply related to the z^* exponent by the hyperscaling relation [ITM 



7 = I - . (6.19) 

This relation turns out to be still valid when D = 1 for closed "1-dimensional membranes" 
with no boundary (self-avoiding loop), but is violated for open "1-dimensional membranes" 
(free open self-avoiding walks). It is also not valid for the physical case of 2-dimensional 
membranes. 

Let us first recall standard results for the partition function of the free membrane 
(6 = 0) (see for instance [^). In this case, we simply deal with the free field, and the 



partition function is given by 



Z(0) = y"l)[r(x)] 5^(ren,)e-^=i['^/^«^ 



det^(-A) ^ 
. Vol(M) . 



(6.20) 



The det'(— A) is the product over the non-zero eigenvalues of the Laplacian on M.^ and 
is defined by analytic continuation in terms of the (^-function C,[s) of the Laplacian as 

det'(A) = e-^'(°) , as) = (6.21) 

n>0 

If one rescales the manifold by a factor X (i.e. the metric by a factor X'^), the volume is 
changed by a factor X^ , and the C function is changed by ({s) C{s)X'^^, and therefore 

det'(A) ^ det'(A)X-2^(°) (6.22) 

C(0) can be calculated explicitly, and is found to be: 

C(0) = — 1 if D not an integer 

— 1-I-— ifD = l, with X = for closed loop, x = 1 for open chain (5.23) 

X 

— 1-1 — ifD = 2 with X the Euler characteristics of Ai 

6 
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The result for D = 2 is valid for smooth open or closed manifolds A4. If Ai has conical 
singularities or if its boundary dAi has corners the result is modified, and can be computed 
using conformal invariance. Thus for the free membrane we get the configuration exponent 

•Jo — 1 = — ud if D not integer 

- ud+ if D = 1 (6.24) 

- ud + -d if D = 2 

6 

with u = (2 — D)/2 the scaling dimension of the free field r. This result has a simple 
interpretation: the l.h.s of ( |6.24| ) is the dimension of the partition function Z{0), —ud is 
the scaling dimension of the constraint 5^ (rem) which removes the translational zero mode; 
the last term, present for integer D, is the contribution of the measure T>[r\ given by the 
trace anomaly. 

We now consider the interacting theory (6 ^ 0). It is sufficient to determine the 
scaling behavior of the normalized partition function 

m - 11 (6.25) 

which can be computed in perturbation theory, using the perturbative rules defined above, 
with the massless propagator GMi^^y) (|6.3|). The constraint rem = fixes the arbitrary 
additive constant in the definition of the propagator so that d^x^y g{x)GM{x, y) = 0. 

We first consider the case D ^ 1,2. Since the renormalized Hamiltonian 7i^[ri^] ( |6.12| ) 
makes the theory UV finite, the renormalized partition function, defined as 

^r(6r;^) = ^ / ^P-(^)] 5^(?Rc™)e-^'^['''^]/'^-^ (6.26) 

is UV finite at e = as a series in the renormalized coupling constant fe^. ii is the renor- 
malization momentum scale. Performing the change of variables ( |5.5|) from renormalized 
to bare quantities, Z-ii^ihn,] fi) is simply re expressed in term of the bare partition function 
2{b) (treating the measure "Dfr] as dimensionless, which is justified in perturbation theory) 

Z{b) = Z{h^)-'''^Z^{h^-ij) (6.27) 

where Z{hii) is the wave-function renormalization factor, and 6(6^, /u) the bare coupling 
constant. If we now rescale the membrane by a factor X, simple dimensional analysis 
shows that 

Z{h-X) = Z{b^{X))-''/^ Z^{b^{X);X^) (6.28) 
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with bji,{X) given hy b = b{bji,{X), X From the renormahzation group flows (|5.7|) and 
we know that for b and n flxed, 6r(X) and obey 



X — Or = — U-— 

dX ^dij 



b^ = ; X^lnZ = 



InZ = 2{iy{b^)-iy) 

^ (6.29) 
For e > 0, 6r(X) — > 6^* as X ^ oo and Z ~ X^*^^ Combining aU these results we see 
that the partition function scales as 

Z{b,X) oc X"^*'^ (6.30) 

hence ( |09|) . 

For D = 1, the same argument applies to the closed chain. For the open chain, we 
must take into account the additional additive counterterm ( |6.15| ). The result is 

7-1 = -z^*rf + 2 77*, (6.31) 

with ?7* the anomalous end-point dimension 

Zq . (6.32) 



d 



b 



For D = 2, unfortunately, the perturbative scheme developed here is useless. Indeed, 
the line D = 2 corresponds to e = 4 and is always at a flnite distance from the line e = 0. 
It is natural to expect that the curvature operator ( p. 13] ) will become relevant, since its 



scaling dimension is always 2 — D, but its mixing with the other operators cannot be 
studied perturbatively by an e-expansion. 

6. 6. Relation with Direct Renormalization 

As a second application of this renormalization formalism, let us derive simply the 
validity of the direct renormalization method used by previous authors. 



We first recall the principle of direct renormalization, along the lines of e.g. [39]. It 
relies on the following ideas: 

(1) One considers only finite membranes, with internal extent L < oo, at e > 0. In 
this case, as we have seen before, physical observables, expressed as perturbative series in 
the "bare" coupling constant 6, are UV and I.R. finite. Then one introduces a physical 
coupling constant g = g(b, L), which is dimensionless and depends on b and on the size of 
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the membrane L. An example of such a couphng constant is provided by the second virial 
coefficient, defined as follows. 

Let us extend the SAM model to the case of two interacting copies and M.' of the 
same membrane. The Hamiltonian is now 



H[r,r']/kBT = \ j {Vvf + \ [ {Vr'f 

^ J M ^ J M' 

b 



+ 



5^(r-r)+ / / 5^(r'-r')+2/ / 5^(r-r' 
M JM Jm' J M' Jm Jm' 



(6.33) 

r and r' describe the embedding of M. and M.' respectively. The partition function Z[h) 
for a single membrane is defined by (|6.17|) , 



Z{h) - j V[y{x)]5\y^^)q-^^^^/^^^ (6.34) 
while the connected partition function for two membranes ^2,c is 

Z2,e(6) = j V[v\V[v']5\v,r^) [e-^P''^']/kB^_e-^['^/kB^e-^[?']/kBTj ^ ^gg^^ 

and the radius of gyration R of a single membrane M. is 

The second virial coefficient is then defined by 

g = -^^-'. (6.37) 

It is a dimensionless quantity, and has a perturbative expansion in h of the form 

g{h,L) = cibL' + C2{bL'f + ■■■ (6.38) 

where ci, C2, ... are dimensionless coefficients which do not depend on the size L of the 
membrane, but depend on its shape. 

(2) It then appears that the scaling functions, i.e. the dimensionless observables F(b, L) 
of the model, which can be expressed as functions of the physical coupling constant g only 



F(6,L) = F{g) , 
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(6.39) 



have a perturbative expansion in g which is UV finite when e — 0, contrary to their 
perturbative expansion in h. An example of a scahng function is the effective scahng 
exponent 

v{h,L) = u{g) = i^^ln(R) , (6.40) 

where R is the gyration radius ( |6.36|) . 

(3) In particular, the Wilson function W{g), defined as 



is a regular function of e when e ^ order by order in g. Its expansion reads 



(6.41) 



W{g) ^ eg - a2(e) g'' + Oig"") , < a2(e) = 0(1) as e ^ . (6.42) 

It thus has a zero Q < g* = 0(e) for e > 0, and from the definition of W{g), one sees that 
in the limit of a large membrane, the physical coupling constant g tends towards a finite 
value: 

g{h, L) ^ g* ; L cx) , 6 fixed . (6.43) 

Since from point (2) W{g) has a regular (/-expansion when e — > 0, inverting W{g) = 
gives g* as a finite series expansion in e. Similarly, the scaling functions also have a finite 
large L limit 

F{b,L) = F(^*) ; L ^ oo , (6.44) 

which correspond to scaling exponents for large self-avoiding membranes. In particular, 

= u{g*) = lim jy{b,L) (6.45) 

L — >oo 

corresponds to the scaling exponent u* obtained in Sec. 5.1. Since F{g) has a regular 
expansion in g when e ^ 0, and since g* has a regular expansion in e, all the scaling 
exponents F* have also a finite expansion in e. 

We now show why direct renormalization is valid. The key property that validates 
the direct renormalization method is (2), namely the fact that the scaling functions have 
a regular expansion in the physical coupling constant g when e ^ 0. This is a simple 
consequence of our renormalization formalism applied to finite membranes. Indeed, one 
can show that the same counterterms as those which make the theory of a single membrane 
UV finite when e ^ also make the model ( |6.33| ) of two interacting membranes UV finite. 
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In other words, re-expressing the bare quantities r and b in terms of renormahzed ones 6^ 
and according to ( |5.5| ) makes the perturbation theory in finite order by order when 
e ^ 0. This property is not completely trivial, but before proving it, let us show that it 
ensures property (2). Therefore let us assume that if we perform the change of variable 

m 

b ^ b^Ztib^){Z{b^)y/^^i' « 6^ = 6^(6^-'^), (6.46) 

the dimensionless scaling functions F have a finite perturbative expansion in fe^ (when 
e ^ 0). By simple dimensional arguments, they depend on the dimensionless coupling 6^, 
on the size of the membrane L and on the additional renormalization scale n only through 
6r and through the dimensionless product iiL. They can therefore be written as 

F{b,L) = F{b^,fiL) , (6.47) 

where F has an UV finite perturbative expansion in b^. Since fi can be chosen arbitrarily, 
we can take /j, = and we obtain that by the change of variable 

6^(6, L) = b^{bL') (6.48) 

the scaling functions can be written as an UV finite perturbative expansion in 6r(6, L) 

F{b,L) = F{b^{b,L),l) = F{b^{b,L)). (6.49) 



This is true in particular for the second virial coefficient (|6.37| ) g{b, L), which is a scaling 
function, and can be written as 

g{b,L) = g{b^{b,L)) = gib^ + 0{bj) ; < r/i = 0(1) when e ^ . (6.50) 

Inverting this relation we can write 6r(6, L) as a series in g, which is UV finite order by 
order when e ^ 0. 

6^(6, L) = bMb,L)) . (6.51) 



Using ( p.49| ) and ( |6.51|) one sees that any scaling fimction F is UV finite when expressed 



in terms of g 

F(6,L) = FibM) = F{g) . (6.52) 
This is nothing but property (2). 
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The Wilson function ( |6.41|) is also a scaling function and is therefore UV finite as 
a function of g. Using ( |6.48| ) , (|6.46| ) and (^) one sees that L-^hn,{h,L) = —I3{hj^), 
and therefore the Wilson function W{g) is simply related to the renormalization group 
/5-function (3{b^) by 

^(9) = - ■ (6.53) 

Thus the renormalization group (7-flow driven by W is just the mirror of the ^R-flow driven 
by /3, and the effective coupling constant g* = ^(6^*), where 6r,* is the I.R. fixed point. 
Similarly, one can show that the scaling exponents given by ( |6.45|) and by (|5.12|) coincide. 



It thus remains to show that the two membranes model ( |6.33| ) is renormalizable by the 
same counterterms ( |5.5|) as the one membrane model. This model is in fact a particular 
case of the general model describing two - a priori different - membranes Aii and A42, 
whose respective configurations are described by the fields ri and r2, interacting by the 
Hamiltonian 



^[ri,r2]/kBT 



1 

2 

2 

62: 
2 



/ d^x, (Vri)^ + I [ 


d''x2 (Vr2)2 


JMi ^ JM 


2 


- [ d^x [ d^y6'^{ri{x 




Jmi Jmi 




- [ d^x I d^y5'^{Y2{x 


)-?2(y)) 


J M2 J M2 




if d^x f d^yd'^{ri{x) 


-Mv)) 


J Ml J M2 





(6.54) 



Indeed, this Hamiltonian is the most general allowed by power counting which describes 
two interacting membranes. Physical observables have a series expansion in the coupling 
constants 611, 622 and 612, which involves the free propagators and Gm2 ^^^J^ since 

the free correlator between the two membranes vanishes 

(ri(x)ri(y))^ = GMi{x,y) (ri(x)r2(y))^ = ^^(^(x,^) {ri{x)r2{y))^ = 

(6.55) 

This perturbative expansion involves the same diagrams as for the one membrane model, 
but the end-points of the dipoles may belong separately to A^i or A^2- It is a not very 
difficult task to repeat the analysis of the short distance singularities of Sec. 3. One finds 
that the UV singular configurations are given by the same molecules as depicted in fig. ^ 
but the points within each atoms, i.e. the subset of collapsing end-points, must belong 
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(1) 



(2) 



(3) 





fiK. 11 



to the same membrane (which depends of course of the atom). One can then repeat 
the MOPE analysis of Subsec. 3.2. Indeed, the coefficients of the MOPE are obtained 
from the short distance OPE of the vertex operators ( 3.13 ) within each atom P, which is 
characteristic of the local properties of each manifold separately. The final result is that 
the short distance behavior of products of multi-local operators ■ ■ ■ ,xp} involving 

the two membranes is given by the MOPE, and that while the operators that appear in the 
MOPE do now depend also on which particular membrane each point Xi of the operators 
belongs to, the coefficients C{y} of the MOPE do not, since all the points in each atom 
belong to the same membrane. The power counting discussed in subsec. 3.3 remains also 
valid, and therefore one deduces immediately that only the molecules with one or two 
atoms are divergent at e = 0. There are now two classes of divergent diagrams with 1 
atom, and 3 classes with 2 atoms, depicted on fig. |TT|. 

These are exactly the same diagrams as those considered in sec. 3. They contribute to 
the perturbative expansion just by different powers of fen, 622 and 612. As a consequence 
the renormalization analysis can also be performed along the lines of sec. 4. Diagrams (1) 



and (2) of fig. |Tl| contribute to a multiplicative renormalization of the first two Gaussian 
terms of the Hamiltonian ( |6.54 ), i.e. to a wave-function renormalization; diagrams (3), (4) 
and (5) contribute to a multiplicative renormalization of 611, 622 and 612 respectively, i.e. 
to coupling constants renormalization. Moreover, simple inspection of the diagrams shows 
that in the perturbative expansion, the divergences (1) and (3), associated only to the 
membrane A^i, depend only on 611, not on 622 and 612. The same is true for divergences 
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(2), (4) and 622- This implies that the general two membrane model ( 6.54 ) is made UV 
finite at e = by the renormalized Hamiltonian 

^ JMi ^ J Ml 

JMi JMi 



+ 
+ 



2 



2 



J M2 J M2 

+ 6r12^6(6rii,6r22,&R12)/"' / X j d^y5'^{Y^i{x) -Y^2{y)) 

J Ml J M2 

(6.56) 



with the same counterterms and Zhip-n) as those of the one membrane theory ( |5.1|) . 

Only the renormalization factor Z{bii, 622? ^12) for gi2 is new. 

A simple inspection of the diagrams shows that if the three coupling constants are 
taken equal b^n = 6^22 = ^R.12 = ^r.: then the three classes of 2-atom diagrams contribute 
with exactly the same factors, and therefore the different Zi, coincide 

Zb{b^) = Zhib^,b^,b^) . (6.57) 

This completes the proof that the two membrane model (|6.33[) is renormalized by the same 



counterterms than the one membrane model. 

Finally, let us remark that a simpler argument (in fact equivalent) is to consider again 
the proof of renormalizability for the one membrane model, and to check that at all the 
steps, the fact that the membrane is connected or not is never important. 

6. 7. Contact Exponents 



As a third application, let us discuss the contact exponents [18] for self-avoiding 



membranes. We refer to for a general presentation of this problem. They are defined 
as follows, as a generalization of the notion of contact exponents for polymers . Let us 
consider the probability that two given points xi and X2 on the membrane are separated 
by the vector r in external space. This probability corresponds to the observable 

P(r;a;i,X2) = ((J'^P - (r(xi) - r(x2)]) . (6.58) 

For an infinite membrane it will have the scaling form 

P(r;xi,X2) = R^2^F{\r\/Ri2) ; i??2 = {[r{xi) - ?{^2)?) , (6-59) 
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where -F is a scaling function. The contact exponents 9 are defined by the smaU r behavior 
of F 

F{r) ~ , r ^ . (6.60) 

For open polymers there are three different contact exponents, depending on the position 
of the points 

9q if xi and X2 are the two end points. Then 9o is given by the scaling relation 6*0 = 
(7-l)/z.. 

9i if xi is a end point and X2 is on the chain 

^2 if Xi and X2 are both on the chain. This exponent is the same for closed polymer 
chains. 

Equivalently, they describe the behavior of P(r; xi, X2) for fixed r, |xi — X2I ^ cxd 

P(r;xi,X2) ~ R^2^-'^ ~ |xi-X2r"(^+^^ (6.61) 

In this limit, we can take r = and hence the scaling exponents 9 are related to the scaling 
dimension of the two point operator 

4>{xi,X2) = 6'^{r{xi) - r{x2)) , {(p{xi,X2)) ~ |xi - X2|^''' |xi-X2|^oo 

(6.62) 

by the relation 

= -u{d + 9) . (6.63) 

The scaling dimension A^ can be calculated in perturbation theory. In our framework, 
perturbation theory is made UV finite when e — by the renormalizations ( |5.5| ). However, 
in order to make the insertions of the bi-local operator (f){xi,X2) UV finite, an additional 
multiplicative renormalization is required, defining the renormalized operator (p^ 

Mxi,X2) = Z46,)5'^(r,(xi)-r,(x2)) = Z^{b^) Z{b^y/^ 6''{r{xi) - r{x2)) . (6.64) 

Z the wave- function renormalization factor and the renormalized coupling constant. 
The scaling dimension A^(6ii) of (p^^ is therefore given (using standard arguments, similar 
to those of subsec. 5.1) by 



d 

A<^(6r) = -dh'{b-n) + /X— ln(Z^ 



= -du{b^) + (3{b^)^\n{Z^{b^)) , (6.65) 
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fig- 12 




fig. 13 



where is the scahng dimension of r^, as defined in ( |5.8|) . The second term in the 

r.h.s of (|6.65| ) comes from the additional renormahzation Z^. The scahng exponent is 
obtained by taking (|6.65|) at the I.R. attractive fixed point 6r,* 



(6.66) 



In order to compute A^, we must determine Z^j,, that is the short-distance singularities 
that are associated with the insertion of ^(xi,X2). These singularities are are obtained 
from the MOPE, and using power-counting one can check that the divergent diagrams are 
those depicted on fig. |T^ , where interaction dipoles collapse toward the inserted dipole 
(xi, X2). These diagrams are very similar to the two atoms diagrams of fig. § and fig. ^ 
but the divergence depends whether the two points, only one point, or none of the points 
are on the boundary of the membrane. We now give a short derivation of the value for 
these contact exponents at one loop. 

Bulk contact exponent 02'- 

This is the simplest situation. At first order on perturbation theory the only divergence 
associated with the insertion of (j){xi,X2) is when the two end-points of a dipole tend 



respectively toward xi and X2, as described on fig. 13 
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Therefore the singularity wiU be proportional to the MOPE coefficient for the con- 
traction (f)(f) (f), calculated in sec. 5.2. Taking into account the symmetry factors, we find 
that the counterterm Zs is at one loop 



1 + -J (47r) 



d , 

2 2 



/ 27r# A 
\T{D/2)) 



Using (|6.63| ) and ( |6.65| ) we get that at one loop 



2+1 r 



D 
2-D 



2D 
2-D 



(6.67) 



-1+1 (. 



27rf 



2+1 r 



D 
2-D 



\V{D/2) 



2D 
2-D 



(6.68) 



and hence 



{2-D) 



l + (2-D)- 



2D 
2-D 



D 
2-D 



+ 0{e^ 



(6.69) 



This agrees with the one-loop result of 



62 can also be easily related to the two-membranes model discussed in sec. 6.7. Indeed, 



the divergent diagrams for the insertion of ^(xi, X2), as described in fig. |T2|, are exactly the 



diagrams (4) in fig. |Tl], with one of the lines between the two different membranes singled 
out. As a consequence, the renormalization factor Zcp{bjO is related to the counterterms 
for the two-membranes model by 



d 



(6.70) 



The derivative with respect to the inter-membrane coupling b^' allows to single out one of 
the links between the two atoms in the diagram. 

Edge-bulk contact exponent 9i: 

In one of the points (let us choose xi) is on the boundary of an open membrane, the 
dimension A,^ of the operator (p is modified. Since we are interested only in local quantities, 
and since for D < 2 no additional boundary operators play a role, we may consider that 
locally around xi, the membrane may be approximated by a half D-dimensional plane 
n = {x± > 0,xu G R^"^). We have seen that for a free membrane, we must choose 
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fig. 14 

Neumann boundary conditions for the field r. Hence, the propagator on 7i is not the 
massless propagator G{x, y) in IR^, given by (|2.13| ), but the Neumann propagator 



G^{x,y) = G{x,y) + G{x,y) 



y = y — 2y± is the mirror of y . (6.71) 



The relevant singularity is now associated to the diagram of fig. |I4. 

The corresponding MOPE coefficient is a simple modification of ( |5.22| ) 

cAJx^-y^,X2-y2) = (An)-''/^ [H^{x^,y^) - G{x2 - y2)r''^^ , 



(6.72) 



where 



H^{xi,yi) = ^{G]^{xi,xi) + GN(yi,yi) - 2G'n(xi, yi)) 



1 



2 \xi - yi 



2-D 



1 



2-D 



(6.73) 



{2-D)Sd V"'""^ 2 
since xi is on the boundary. The singular part of such a diagram is similar to that of ( |5.30| ) 
and, denoting u = yi — xi and f = 1/2 — 2:25 it is given by the small u and v singularity of 
the integral 



T{D/2) 



'u±>Q 



2(2-D)7r^/2 
4 u^ + ul 



-d/2 



X 



2-n 

— \2u±\ +2 \v\ 



1 -d/2 



(6.74) 



This singular part can be calculated explicitly. It is equal to the singular part ( |5.30| ) (which 
gives the counterterms Zi, and Z^j, in the bulk- bulk case), multiplied by the following factor. 
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fig- 15 

a;(D), which is simply given by the shift from the propagator G to the Neumann propagator 



u{D) 



2-D 



1 



J d'^u 5 - 1] 



s 



D-1 



2S 



D Jo 



du± 



11. 

2 + J 



2-D 



D 
2-D 



11. 

2 + l'^"^' 



2-D 



2 

2-D 



U 



D-3 
2 



with M™' 



[2(1-2- 



D\-\ 2^ 



(6.75) 

and Sd = r{D/2) volume of the unit sphere in IR^. This integral is convergent for 

1 < D < 2, and can be estimated numerically. When D —* 1 uj{D) — ^ 1/2, and it can be 
calculated with a finite part prescription at u± = u™^^ for < D < 1. 

Repeating the argument for 6*2, we find that at order e, the edge-bulk contact exponent 
6*1 is simply 

6*1 = u{D)e2 + 0(e2) (6.76) 



Edge-edge contact exponent 6*0: 

A similar reasoning allows to compute the edge-edge contact exponent ^o- Now both 



points xi and X2 are on the boundary, and the relevant singularity is given by fig. |T5 
with MOPE coefficient 



CNNL(xi-|/i,X2-y2) = (47r)-^/2[ifN(xi,2/i)+if: 



N[X2,y2j 



-d/2 



(6.77) 



The final result is at first order in e 



6*0 = u{D)''e2 + C(e2) 
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(6.78) 



Comparison with previous results: 

These resuhs for 9o and 9i differ from those of |TH] for an open membrane with 



a smooth boundary. Indeed in one obtains at order 0{e) the result 62 = 26 1 = 
46*0, with 62 given by (|6.69|) , which differs from (|6.76| ) and (|6.78| ), except for the case of 
open polymers D = 1 (since uj{1) = 1/2). This discrepancy in easily explained by the 
difference in the choice of boundary conditions. In [|T^ , an infinite membrane A4 ~ IR^ 
is considered, with self-avoidance interaction acting only for points within a finite domain 
T> G Ai (with a smooth boundary). For points outside T> there is no self-avoidance, so 
the outside membrane is phantom-like. This boundary condition turns out to lead to 
different edge exponents than our physical Neumann boundary conditions, which describe 
a finite membrane Ai (with a smooth boundary), with self-avoidance interactions. This 
fact should not be considered as surprising. 



7. Other Models 

7.1. Long Range Interactions 

The MOPE that we have derived for the multilocal operators ^{A,m}{2;} given by 
( |5.10| ), involving the contact interactions d'^{r{x) — r{y)), can be generalized to other op- 
erators, such as the long range Coulomb interaction (in external d-dimensional space) 

^c{x,y) (X |r(x)-r(|/)|2-^ . (7.1) 

It has a representation similar to ( |2.8| ) 

^c{x,y) = I^!^^^d\k^+k2)\k^\-'e'^^-'(^^e^^^-'^y^ (7.2) 

The only difference with the contact interaction is the appearance of a new dipolar con- 
straint replacing (|2.9|) : 

Cc{ki,k2} = (27r)'^5'^(ki+k2)|ki|-2 . (7.3) 

The short distance behavior of a product of such Coulomb operators can be analyzed 
according to the analysis of Subsec. 3.2. The corresponding MOPE will generate the same 
contact operators ^{A,m}{x} as in ( |3.25| ), with different coefficients C^'^''"-^. Indeed, the 
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original dipolar constraints C, now replaced by Cc, appear only in the coefficients ( |3.22|) , 



not in the operators (|3.23 



If we now consider charged tethered membranes (with internal dimension D) with 
non-screened repulsive Coulomb interactions in d dimensional space, which generalize poly- 
electrolytes, they are described by the Hamiltonian ( |2.1| ), with d'^{r{x) — r(x')) replaced 
by "i^cix, x'), and a new coupling constant 6 > corresponding to the squared charge den- 
sity. The consequences of the new MOPE are twofold. (1) The power counting is slightly 
modified. The Coulomb interaction is relevant at large distance if 

e' = 2D - {d - 2)^--^ > 0. (7.4) 

For e' close to 0, only the Gaussian operator (Vr)^ and the Coulomb interaction need 
to be renormalized. (2) Since it is not generated by the MOPE, the long range Coulomb 
interaction is in fact not renormalized. This result holds to all orders in perturbation 
theory. Physically it means that there are no screening effects for a uniformly charged 
fluctuating object. 

As a consequence the renormalized Coulomb Hamiltonian is of the form 

^^AbT = ^ I d^x (V.r,(x))' + ^ I d^x' |r,(x) - r.(x')|^-^ , 

(7.5) 

without coupling constant renormalization factor Z^- This renormalized Hamiltonian is 
equal to the bare one by the change from renormalized to bare quantities as in (|5.5| ) 

r = Z^l\h^)v^ ■ b = 6,Z(^-2)/2(6^)/ . (7.6) 

Therefore the RG function /9(6r) and the anomalous dimension z^(6r), given by ( |^ ) and 
( p78|) , are not independent, but related by 

/3(6r) = -2Db^ + {d-2)b^u{b^) (7.7) 

If, as can be checked at one loop, an UV fixed point 6^* > exists, such that /3(6r*) = 0, 
then the anomalous size exponent u* = u^b^*) is exactly given by 

= ^ . (7.8) 

d-2 ^ ' 

This result coincides with the result obtained using a Gaussian variational method 
40| - |4^ , an approximation which in this case turns out to be exact. This result is of course 



meaningful only if u* < 1. If the r.h.s of ( [7.8|) is larger than 1, this means that the long 
range Coulomb interaction makes the object rigid and completely flat, with u = 1. This 
occurs when d < 4 for polyelectrolytes {D = 1) and when d < 6 for charged membranes 
(D = 2). 
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1.2. Tricritical Behavior: versus G" Points 

Another application concerns interacting tethered membranes at the tricritical point. 
One might expect on physical grounds that for "real" tethered membranes, there is a 
competition between two-body attractive interactions (for instance long-range Van der 
Walls forces) and hard core repulsive interactions (for instance hydration forces). This is 
known to occur for polymers, and as a consequence there is a low temperature collapsed 
phase, where attractive forces dominate, and the high temperature swollen phase, where 
repulsive forces dominate. These two phases are separated by a critical O point. In the 
framework of the Edwards model, the effective 2-body coupling b vanishes at this point, 
and one must introduce the next most relevant operator. For polymers it is well known 
that this operator is simply the 3-body repulsive contact interaction [53|T^]. 

For membranes with arbitrary 1 < D < 2, the situation is more subtle. Among all 
multifocal operators ( |6.1UD one has to consider two potentially relevant operators: the 
3-body operator 

$3{x, 2/, z} = d^{r{x) - r{y)) d''{r{y) - r{z)) , (7.9) 
with engineering dimension £3 for the corresponding coupling 

63 = 3D-2d^^^^ ; (7.10) 
and the modified 2-body operator 

^'^{x,y} = -Ar-5''(r{x)-r{y)) (7.11) 
where is the Laplacian in (i-dimensional space, with corresponding dimension €2 

e'^ = 2D -id + 2)^^^ . (7.12) 

This last operator describes a 2-body interaction repulsive at short range and attractive 
at larger range. 

A simple analysis shows that the critical ©-line 63 = and the 0"-line €2 = intersect 
at (D = 4/3, (i = 6). If D < 4/3, and in particular for polymers, the 3-body operator j fl^ ) 
is the most relevant, while for D > 4/3 the modified 2-body operator ( |7.11| ) dominates. 
Considering only this last operator, the one loop calculation of the renormalization group 
fiow and of the scaling exponent u* at order 62 is presented in the next subsection. The 
corresponding one loop results for the 3-body operator require numerical calculations, 
and the analysis of the cross-over between the two regimes at L> = 4/3, (i = 6 is more 



complicated to study. They can be found in [44 
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7.3. 1-Loop Results for the 6" Point 

We first apply the MOPE formalism to compute at one loop the renormalization group 
functions for the O" point. We need the MOPE for $2(^15 ^2} local operators when 
xi X2, which is 

<5'^{xi,X2} = Cl„{xi-X2}l{x) + Cl?,{xi-X2}:Vc.r{x)V(3r{x): + ... (7.13) 

with X = (xi +X2)/2. Starting from the OPE ( p.l6| ) and integrating over the k, now with 
an extra kernel k^ to account for the — Ap, we obtain 

C%|{x^-X2} = -^(^ + l)(47r)-^/2(-G(xi-X2))-2-^/2(x«-x«)(xf-x?) (7.14) 
The MOPE for $'2'*^' ^ is of the form 

$^'{xi,yi}$^'{x2,2/2}^ =^ Cj?!^„{xi -X2,yi -y2}$2{2;,2/} + ■■■ (7-15) 

Ki^y2 

where x = (xi + X2)/2, y = {yi + 2/2)/2. We use the OPE ( |5.2C1| ) and obtain, along the 
same line as in sec. 5.2 



C|?,^,,{xi-X2, 2/1-2/2} = i^n)-"/' + [_G'(xi-X2)-G(2/i-2/2)]-'-'/' 

(7.16) 

We are now in position to compute the renormalization factors. The renormalized 
action is 



2 



At one loop, one obtains Z" from the singular part of the integral of C^^{xi — X2} 



as e'4 0. 



Similarly we obtain Z'^ from Q>^i^n 



2 ^2 

2 



e^' 2 



^6+8 2;l^r(D/2); r(^) ^ ^ 
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The renormalization group functions (3 and u are computed using 
I. We find for the anomalous exponent at the G" point (at order £3) 



(|5.7|) and 



2-D 



V 



+ e 



4D 



D 
2-D 



{2-D) 



2D 
2-D 



dj_ I 3d _ 1 
16 8 2 

(f + 1) 



n -1 



(7.20) 



8. Conclusions 

In this paper we have given a general formalism to study in perturbation theory the 
generalized Edwards model, which describes the continuum limit for fiexible self-avoiding 
tethered membranes and polymers. We have shown that the short distance behavior of this 
model in encoded into a multilocal operator product expansion (MOPE), which generalizes 
Wilson's operator product expansion. This MOPE allows us to study in a systematic way 
the short distance behavior of the model, and to prove that it is renormalizable at all orders 
in perturbation theory. The formalism constructed here allows in particular to prove the 
validity of the direct renormalization method, which was used previously to derive at first 
order in perturbation scaling laws and scaling exponents for large self-avoiding tethered 
membranes. It also provides a direct proof of the consistency of direct renormalization 
for polymers. We have also presented a few applications of this formalism, including 
membranes with long range interactions and membranes at the 0-point. 

Let us stress that our formalism is quite general, and should apply to a large variety of 
interactions, provided that they are multi-local (with respect to the "internal" membrane 
space), singular at short distances in the external embedding space, and invariant under 
global translations. The method applies also naturally to the case of several interacting 
membranes and polymers, and might be applicable also to membranes in disordered media, 
via the replica method (as already done for polymers). 

Let us also stress that the method is essentially perturbative. Several important 
issues still have to be addressed in this context. Is the perturbation series asymptotic, 
and in which sense is it summable? In Quantum Field Theory it is expected, for instance, 
that the e = 4 — d expansion is Borel summable; this conjecture appears to be satisfied 
numerically, and is the basis of very efficient ressummation methods which lead to very 
accurate predictions for critical exponents. Almost nothing is known, even at the heuristic 
level, for the generalized Edwards model. This is a very important point, since to reach 
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the point of physical interest D = 2 one must extrapolate up to e = 4. In this respect, it 
would be suitable to control the large orders of the perturbative expansion. 

One thing which is known is that at the first and second orders, the estimate for 
the size exponent u coincides at large d with the result z/var = 2D/d of the Gaussian 



variational method |j2^,^. Arguments have been given that this is true to all orders in a 
1/d expansion ||2^, with a difference z/ — u^ar = 0(exp(— (i)). One should notice however 
that this property seems to be true only for the size exponent u. Also one should stress 
that no systematic 1/d expansion has been constructed so far for this model. 

Another important issue is the relevance of curvature terms. It is known for two 
dimensional objects both the intrinsic (or Gaussian) curvature and the extrinsic (or mean) 
curvature might become relevant. For polymerized membranes, the intrinsic curvature does 
not fluctuate (see the discussion of section 6). So it should not play an important role, 
except for the configuration exponent 7, which is expected to get an extra dependence 
on the integrated curvature, i.e. of the Euler characteristics, for D = 2 exactly. On 
the contrary, the size exponent u is insensitive to the intrinsic curvature. More drastic 
effects are expected from the extrinsic curvature, i.e. from the bending rigidity of the 
membrane. For phantom membranes (without self-avoidance) but with bending rigidity it 
is known that a crumpling transition occurs, separating a flat phase (at high rigidity) from 
a crumpled phase (at low rigidity). In the flat phase self-avoidance is clearly irrelevant. The 
generalized Edwards model considered in this paper is applicable to the crumpled phase, 
and concentrates on the modification of its properties due to self-avoidance. It implicitly 
assumes that bending rigidity terms are irrelevant, which is true close to e = (by power 
counting) . The important point is their relevance at the infra-red stable fixed point (which 
describes the crumpled phase) for finite e, in particular for e = 4 {D = 2). A heuristic 
analysis can be obtained by simply comparing the size exponent u in the crumpled phase 
(with self-avoidance) with the size exponent t'cr of phantom membranes at the crumpling 
transition. It is expected that bending rigidity is relevant whenever 1/ > Ucr, and that this 
will drive the membrane into the fiat phase The different estimates of and 



of lead to a lower critical dimension di^ below which two dimensional self-avoiding 
membranes should be fiat, lying in the range 3 < die < 4. Clearly it is a major issue to 
have more refined estimates of d\c, in order to establish whether physical two dimensional 
membranes in three dimensional space may be in a crumpled state or are always fiat, a 
behavior which would be very different from that of ordinary polymer chains. 



85 




Acknowledgments: 



fiK. 16 



F.D. thanks K. Wiese for his interest and useful remarks. We thank M. Bergere for a 
careful reading of the manuscript. 



Appendix A. General form of the MOPE 

It is not difficult to generalize the multilocal operator product expansion (|3.25|) , which 
describes the short distance behavior of a product of bilocal operators 5'^(r(x) — r(y), to 
describe the short distance behavior of a product of the most general multilocal operators. 
These operators $ are given by ( |3.24| ), they are characterized by a set X of "end-points" 
by local operators Ai and by (i-uples liii G IN'^: 

$ = = y'd'^r n{:A,(xO(iVr)'^^5'^(r-r(a:0):} . (A.l) 



We depict graphically such an operator by the starfish-like diagram of fig. [16. 

We now consider P such operators $i, ■ ■ ■ $p, with different end-points, and we are 
interested in the behavior of the product of these operators when we partition the set of 
all end-points X = Xi U ■ • • U Xp into a set of atoms V , assign a given position x-p to each 
atom P, and let the end-points Xi, i & V inside each atom V tend toward x-p. In this limit, 
the operators form a molecule M., that we assume to be connected. 

The general MOPE takes the form of a sum over all multilocal operators ${a;p} 
relative to the set of atoms V: 

$i{xa---$p{xj = 5^C|^...^^{yJ-${a;p} (A.2) 
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where the coefficients C| ...^ depend on the relative coordinates {yi }ie'P inside 

each atom V, and not on the coordinates x-p of the atoms themselves. They are homoge- 
neous functions of the y^'s, with degree degC which depend on the operators $i, ■ ■ ■ $p, 
of $ and of the molecule Ad. 

The derivation of (|A.2| ) is parallel to that of ( |3.25|) . One starts by rewriting each $j 



/ = 1, ■ ■ ■ , P as 

J J ieii (27r)'» ieii 
The integration over rj gives a neutrality constraint for each starfish 

Ci{k,} = {2nf6''{ E k,) (A.4) 
ieii 

Inside each atom P, the product of vertex operators dressed by the AiS has an operator 
product expansion (which generalizes ( |3.13|)(|3.14|)(|3.15|) ) of the form 

n ■.Mx,)e^^^'^-^^: = VC7r\U2/„k,}e -^^^ :A(a:p)e : (A.5) 

ieJi A 

where the coefficients C^^ y are monomials in the k^'s, and depend only algebraically 
of the y^'s, i.e. contain only (possibly non-integer and negative) powers of the rela- 
tive coordinates i/i — yj. Inserting a momentum h-p for each atom via the identity 
1 = / (i'^k-p5'^((Eiep ki) ~ kp), using the connexity of the molecule Ai to rewrite one 
of the starfish neutrality constraints X/ as a global neutrality constraint for the molecule 
5'^(E-pk-p), and now Taylor expanding in the k-p one ends with (|A.2| ), with the explicit 
expression for the coefficient 

cl...Ay^}^ /.n|;;^,n'c.{ka n/d-^'x 

(A.6) 



, , 1 . . ^ -5 S k,-kjG{y,,yj] 

veM i-^»/'6T' m-p! ^ iev 



with n' being the product over the /'s but one. 
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Appendix B. The (massless) propagator on the sphere So- 

The propagator for a free field on the D-dimensional sphere Sd with constant curva- 
ture can be defined in various ways which allow for a proper analytic continuation in D. 
Let us consider the unit sphere with radius L = 1 and volume 

D + l 

yo\{SD) = So+i = (B.l) 

The massive propagator G(x,y;t), where t is the squared mass t = m? , is the solution of 
the equation 

{-/\^ + t)G{x,y-t) = — L=5(x,y) (B.2) 

where A^; is the covariant Laplacian. By 0{D) invariance it depends only of the geodesic 
distance I between x and y, G = G{i; t), and in radial coordinates ( p.2| ) becomes 



G{i;t) = for £ ^ mod(27r) (B.3) 



supplemented by the condition that at the origin the singular part of G{i; t) behaves as 
the massless propagator in plat space 



,2-D r(D/2) 



By the change of variable cos(£) = z (|B.3| ) reduces to a Gegenbauer equation 

(1 - z^)G" - DzG' ~ tG = (B.5) 

and G{i; t) is expressible in terms of known special functions (Gegenbauer functions or 
hypergeometric functions). 

In the massless limit t — the propagator becomes I.R. singular. Subtracting the 
I.R. singular part due to the zero mode of A, we obtain the massless propagator 

GM = lim(G(£;t) (B.6) 



which is the solution of (|6.3|) on Sd- 

Explicit integral representations for the propagators on the unit sphere can be obtained 
for instance by the following argument. Let us embed the unit sphere Sd in IR^"*"^ and 
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use the spherical coordinates X — > (r, (j)) = (r, (j), ip\, ■ ■ -(^0-2) in IR'^"'""^. If G{i] t) is the 
massive propagator on iSd, it is easy to see that the function on IR^"*"^ 

F{X-t) = Hir;t)Gi<P;t) (B.7) 

satisfies the equation 

-A^F = H{r;t)r-^S{$) (B.8) 

where d{(p) is the delta function at the origin on the unit sphere, provided that H{r;t) 
satisfies 



^ D d ^ t' 



H{r;t) = . (B.9) 
Hence, from the knowledge of the massless propagator on IR^"*"^, 



\X -Y 



l-D 



we can obtain the function F{X) as the Coulomb potential created by a linear distribution 
of charge p(r) = r^~^H{r; t) along the positive unit axis Y{r) = (r, • ■ ■ 0) 

/"OO 

F{X) = I drH{r;t)r^-^G{X,Y{r)) (B.ll) 

Solving ( [B.9|) leads to 







ii-(r;t) = r — ^y\ — ) "* (B.12) 

and from F{X] t) we obtain finally the integral representation for the massive propagator 
on the sphere 

G{c^-t) = r - r'^Wi^Y-' [r2 + 1 - 2rcos((/>)] 1-^/2 (B.13) 

4 7r 2 Jo ^ 

The lb representations are equivalent via the change of variable x 1/x. The integral 
( [B.13| ) is convergent for 1 < D, and allows for an explicit analytic continuation in D. For 
< 1 it is defined by a finite part prescription. 
The massless propagator G{(j)) is defined by (|B.6|) , and has the integral representation. 



r(^) f°-dr 



D + 1 

4 7r 2 Jo ^ 



l-D 



G{cl)) = '-^-i^ / - [r2 + l-2rcos((/>)]^-^(l-r)) , (B.14) 
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with 9 the Heaviside step function. If we shift G it by a constant so that G{0) = 0, we get 



G{cP) ^ G{cp) - G{0) 



4 TT 2 



^ (^[^2 + 1 _ 2r cos((/>)]^ - [(r - 1)']^) 



FinaUy, from ( [B.13| ) one can construct the heat kernel K{(p; s) on the D-dimensional 
sphere through the inverse Laplace transform with respect to t 



K{4>-s) = / —e'^GiM (B.16) 
^ 2i7r 



Appendix C. The short distance OPE in curved space 

In this appendix we give the outline of the derivation of the OPE for vertex operators 
( |6.8|) in curved space. It relies on the short distance behavior of the free massless propagator 
G{x,x') = (r(x)r(x'))o on a general curved manifold M., when x,x' — z. This behavior 
can be inferred from the small proper-time behavior of the heat-kernel 

K{x,x'-s) = (x|e"^|x') (C.l) 

using the relation 



/"OO 

G{x, x') = ds 
Jo 



Vol(7W) 



(C.2) 



The last term in ( |C.2|) is required in order to subtract the contribution of the zero mode 



of the Laplacian A = D'^Da, which contributes to the large s limit of the heat-kernel: 
lims^ooK{x,x';s) = l/Vol(A^). 

The small s behavior of K is known to be given by the following asymptotic expansion 
(see in particular the chapter 9 of 



oo 



K{x,x';s) = (47rs)-^/2e-^'("'"')/4" ^ a,(a:,x')s^ . (C.3) 

r=0 

£{x,x') is the geodesic distance between x and x' . D is the dimension of the manifold A4. 
A crucial property of the functions i'^{x,x') and of the coefficients ar{x,x') is that, for a 
smooth metric, they are smooth functions of x and x', at least as long as these points are 
close enough. Their Taylor expansion at x,x' = z depends only on the local properties 
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of the metric at z. This means that at a given order m, the coefficients of their Taylor 
expansion are a differential polynomial in the metric, that is involve only the metric tensor 
g{z) and its derivatives dg{z), ddg{z), etc. . . at ^ up to some order (which depends on m). 
For instance, the expansion of the coefficients ar{x,x') starts as 

ao{z,z) = 1 ; ai{z,z) = ^ R{z) 

1 11 1 (^-4) 

a2iz,z) = -D^D^R{z) + ^R^i^) - —R<-0 R^p + _ i?^^^, i?^^^, 

which means that the term of degree of the Taylor expansion of a^- involves exactly 2r 
derivatives d of the metric. Then, the terms of degree m of the Taylor expansion of 
(that is the sum of the terms of order [x — z)^{x' — z)^ , with m = |m| + are 
differential polynomials of the metric which contains exactly 2r + m derivatives d (that is 
polynomial in g~^, g and its derivatives dg, ddg, dddg, etc. . . at the point z). Similarly 
for the squared geodesic distance we have 

i'^(XjX') = (x — x')" (x — x')(3 gaffiz) + ■■■ 

- II - ... ^^-^^ 

1 1 1 1 Z ^ 

and the terms of degree m > 2 contain exactly m — 2 derivatives. 

The (x, x') — > z behavior of G{x, x') is obtained by using (|C.2| ) and ( |C.3|) . For s ~ ||x — 
x'Wz in the integral (|C.2|) the small s behavior of the heat-kernel dominates. This generates 
a singular contribution for the small distance expansion of G, which is proportional to 
||a; — x'\\1~^ , times powers of {x — z) and {x' — z). 



r=0 

J\\2-D 



(C.6) 



I _ inZ-D f-^f l\ 



where C{x, x') is an asymptotic series in x — z and x' — z. The coefficients of order m of B 
are differential polynomials of the metric g at z, with exactly m derivatives d. Therefore 
( |(J.6D can be written as a sum over differential monomials B[g](z) of the metric at z, of 
the form 

G'sing(x,x') = ^C^(x-;2,x'-z))S[(7](z) (C.7) 

B[9] 

where the coefficients are homogeneous function of (x — z) and {x' — z) of degree 
2 — D + m, where m is the total number of derivatives d in B[g]. The leading contribution 
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is given by the r = term, and is associated to the unity operator Sq = 1, it is nothing 
but the flat space propagator Gq{x,x') = — ||x — x'|p~^/(2 — D)Sd at z. 

The rest of the s-integral in ( |C2D generates a regular contribution G-^e^i^x, x'), which 
is an asymptotic series in integer powers oi x — z and x' — z. By a formal Taylor series 
expansion the coefficients of this series define the expectation values of the normal ordered 
operators quadratic in r 



aeg(x,x') = E J^'" (Ar[V"^rV"^'r](z))o (C. 



The total asymptotic expansion of the propagator is 

G{x,x') = Gj-egix,x') + Gsingix^x') (C.9) 

This decomposition of the short distance expansion of the propagator in terms of the non- 
analytic part Gsing and the analytic part Greg is unambiguously defined if D is not an even 
integer. This short distance expansion is nothing but the Operator Product Expansion for 
two free field operators in curved space. The operators that are generated are the normal 
ordered operators AAfV^rV" r] and the purely geometrical local operators B[g]. 



Now we can derive the general OPE ( |6.8|) . The product of exponential operators can 



be separated as in (|3.13| ) into its normal product : [ ] : times its expectation value 

■|-|-gik,.r(x,) _ ^-Qgik,.r(x,)^^. J-|-gik.-?(xO. (C.IO) 

i i i 

The normal product : [ ] : is defined in the usual way by subtracting all propagators between 
pairs of x^'s, and it is therefore analytic in the x^'s when they tend toward the same point 
z. The singularities when Xi ^ z are contained into 

_^ "2 ^ ^ -^i ' Gi^Xi ,X j ) "2 ^ ^ ' Greg (^i i ) "2 ^ ^ ' Casing (^i 1 ) 

(]^e''^^'''^^^^)o = e = e e 



(C.ll) 



Let us define the normal ordered product 



n 



giki-r(xi) 



2 ^ ^ "■^J Greg, (^i ■i^j ) 
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Since Gsing{xi, Xj) has an asymptotic Taylor expansion, we can Taylor expand this normal 
ordered product around z in powers of the Hi = Xi — z. We thus obtain an OPE similar to 



iki-r(xi) 



Air] 



(C.13) 



(k = hi) where the sum runs over the local operators A[r](2) which are monomials in 
the derivatives of r, that is on the same operators than in (|3.15| ). Since they are obtained 
from a formal Taylor expansion, the coefficients are also exactly the same than those 
in (|T5D. 

Now we separate the singular part of the propagator into 



Gsingixi.Xj) = Go{xi,Xj) + ^ C^ivi^yj) B[g]{i 

B[g]^Bo 



(C.14) 



and expand the last, subdominant sum in the exponential (|C.11|) . We obtain also an OPE, 
which involves only geometrical local operators B[g] 

2 ^ ^ '^j Casing i^i i^ j ) "2 ^ ^ '^3 ^0 {^i ) 

e - 5^C^{2/„kai?[(7](z)e (C.15) 

Big] 



Putting together the OPE ( |C.13|) and (|C.15| ) we obtain an OPE of the form 

^ kt-kj Go(xi,xj) 



(C.16) 



where now the sum runs over all products A = A[r\B[g], that is over all local operators A 
which are monomials in derivatives of r and depend on the metric g and its derivatives at 
z. The coefficients are nothing but = C^^^C^^^'^ . This ends the derivation of the 
OPE in curved space. 
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Figure Captions 



Fig. 1. 
Fig. 2. 
Fig. 3. 

Fig. 4. 

Fig. 5. 

Fig. 6. 
Fig. 7. 
Fig. 8. 

Fig. 9. 
Fig. 10. 
Fig. 11. 
Fig. 12. 
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Fig. 14. 
Fig. 15. 
Fig. 16. 



A dipol^ 
A Diagram 



A singular contiguration| is obtained by contracting some points and by forming 



"molecules" . 



General structure of UV divergent diagrams with = 1 (a) and N = 2 (Jo) 
atoms. 

Rooted UV divergent diagrams with = 1 and N = 2 connected components 
(atoms). 



The diagrams giving [divergences at one loop 
The diagrams giving [divergences at two loops] 

The 1-loop diagrams containing insertions of 1-loop counterterms and giving also 



divergences at two loops] , with the corresponding nested 2-loop divergent dia- 
grams. The cross denotes the insertion of a local 1 or (Vr)^ operator. 



Tableau construction in the simple case of 2 diagrams 



The diagrams giving [divergences at one loop 

Divergent diagrams for the 2 membrane] model 

Divergent diagrams for the insertion of $(xi, X2) 

L)ne loop divergent diagram for the insertion of X2) 

L)ne loop divergent diagram when one point is on the boundary. 

(Jne loop divergent diagram] when both points are on the boundary. 



Starfish representation of ]a multilocal operator 
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